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WE'RE REACHING INTO SPACE 


Bell Laboratories research 
with chilled ruby amplifiers 


speeds the day we may 
telephone via satellites 


A strange combination of Nature’s forces at Bell Laboratories fore- 
shadows the day when world-wide phone calls may be relayed via man-made 
satellites orbiting the earth. It is a union of synthetic rubies and extreme 
cold, making it possible to amplify microwave signals from these satel- 
lites clearly. 


Synthetic rubies possess an extraordinary property when deeply chilled 
and subjected to a magnetic field. They can be excited to store energy at 
the frequencies of microwave signals. As a signal passes through an excited 
ruby, it releases this energy and is thus amplified a thousandfold. 


Bell Laboratories scientists chose a ruby amplifier because it’s 
uniquely free of “noises” that interfere with radio signals. For example, 
it doesn’t have the hot cathodes or hurtling electrons that generate noise 
in conventional amplifiers. It is so quiet that there remains only the noise 
made by matter itself in heat vibrations. But at a temperature close to 
absolute zero, this also is silenced. Even very faint signals from satellites 
can be clearly amplified and studied for their practical possibilities. 


Bell Laboratories scientists were first to discover that matter itself 
generates electrical noise. They also discovered that stars send radio waves, 
and thus helped found radio astronomy. It is particularly fitting that the 
same scientists, in their endless research on noise, should now battle this 
number-one enemy of telephony in the dramatic new field of communica- 
tion via satellites. The ultimate goal, as always, is the improvement of your 
Bell System communications services. 


BELL TELEPHONE LABORATORIES 


t/ World center of communications research and development 
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ON THE SUMMATION OF > fn AND ITS ASSOCIATED 
n=0 
INTEGRALS * 


BY 
E. I. JURY! AND M. A. PAI! 


SUMMARY 


A geometrical interpretation of the convolution z transform is given and this is 
applied to develop z transforms of functions of the type f(?)/t*, k > 0. After estab- 
lishing the equivalence between certain forms of integrals, the authors have used the 
integrals to sum certain forms of infinite series. Possible extensions of this technique 
to other forms of series and to sum finite series are also indicated. 


INTRODUCTION 


z transformation technique has been successfully applied to solve 
linear difference equations. Convolution of z transform and modified 
z transform developed by Cypkin (1)? and one of the authors (2), re- 
spectively, will be used in this paper to derive the z transform of certain 
discrete time functions of the type (f,)/(n*). These z transforms will 
then be utilized to sum infinite series, which are well known in mathe- 
matical literature. During the discussion equivalence between certain 
types of integrals and their use will be established. First the complex 
convolution of two discrete functions will be explained constituting a 
brief review of (2), and the paths of integration will be defined. Their 
importance will be evident later. 


1. COMPLEX CONVOLUTION (1, 2) 


Suppose we have two continuous time functions f(t) and A(t) having 
F(s) and H(s) as their Laplace transforms. Let f*(¢) and h*(t) be the 


* This research was supported by the U. S. Air Force through the Air Force Office of 
Scientific Research of the Air Research and Development Command, under Contract No. 
AF 18(600)-1521. Reproduction in whole or in part is permitted for any purpose of the 
United States Government. 

1 Electrical Engineering Department, University of California, Berkeley, Calif. 

2 The boldface numbers in parentheses refer to the references appended to this paper. 


(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
the JourNAL.) 
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sampled functions, their Laplace transforms being F*(s) and H*(s). It 
is required to find the Laplace transform of [f*(t)h*(t) ]. 
Using complex convolution in the s plane (2, 3) we have 


e+j(r/T) 


LE f*(Dh*(t)] = F*(q)H*(s — g)dq. (1) 


c—j(r/T) 
Both F*(s) and H*(s) have an infinite number of poles in the s plane 
because of the sampling action. Assuming as an example that F(s) has 
a pole in the left hand plane at s = s,, and H(s) has a pole also in the 


left hand plane at s = ss (both s, and sg real), then the pole-zero pattern 
in both the s and g planes of F*(s) and H*(s) is as shown in Figs. 1, 2.° 


q-PLANE 
-PLANE x Im x 


x 


er 
T x x h | 


am \ Re 
x x 77 \ x 
x x x T x 
T x SB = 

Fic. 1. Pole-zero pattern of F*(s) and H*(s) Fic. 2. Region of analyticity in the q plane 


in s plane. and path of integration. 


The strip of analyticity in the g plane is shown cross-hatched and the 
path of integration is shown dotted. Since both F*(s) and H*(s) are 
functions of e7* only, we can replace e™* by z and e™ by p so that F*(q) 
becomes a function of » only and H*(s — g) becomes a function of 2/p 
only. As required by the mapping dg = (1/(Te™))dp = ((p-)/T)dp. 
The mapping of the s plane into the z plane and g plane into the p plane 
should be obvious by referring to Figs. 3 and 4. Note that the strip of 
analyticity in the g plane maps into an annulus of convergence into the 
p plane. With these substitutions 


where I is the dotted contour in the p plane. 

By Cauchy’s theorem on residues we can evaluate residues at the 
singularities of p-! F*(p). If the value of the integral along the infinite 
circle tends to 0 then we can also evaluate residues at the singularities 


3 Equation 1 is applicable even if F(s) and H(s) have poles in the right-hand plane. 
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(2) of H* (<). Therefore 


1 
>< Residues at singularities of p-'F*(p) 


— > Residues at singularities of H* ( 5 ). 


Z7-PLANE 


—PLANE 
P 
UNIT CIRCLE Nee INFINITE CIRCLE 


UNIT 


Rapius 


(CONTOUR 
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Fic. 3. Poles of F*(z) and H*(z) in z plane. Fic. 4. Path of integration in p plane. 


This can be illustrated by a simple example. 
Consider: 


Zz 
H*(z) 


F*(z) = 
1 
2njJrp—az— 
Jr 


aL f(t)h(t)] = 


Denominator of the integrand is at least two degrees higher than the 
numerator in p. Hence the value of the integral along infinite circle is 
zero. 

Now evaluating residues inside the contour I’, we have 


al = ——. (5) 
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81 
= 
rf 
Tm 
—— = 
= 
Now 
ee 
dp 
pb 
| 
Z 
dp (4) | 
p b 
me 
: 
4 
: 
: 


82 E. I. Jury anp M. A. Pat (J. F. 1. 


Evaluating residues with respect to singularities of H* (2). we have 


(6) 
a 


z—ab 


2. z TRANSFORM OF f(t)/ (k IS AN INTEGER > 0) 


The z transform of f(t) /t has been derived in a previous paper (2) as 


1 n 


provided f, = 0 and lim f,/n exists. We will show that this result is 
n—0 n—0 


equivalent to 4{f(t)-1/t} obtained by using convolution. Using the 
above formula we first of all obtain 4(1/t) fort > 0. If f(¢) is a unit 
step function starting at ¢ > 0 then F*(z) = 1/(¢ — 1) so that applying 
the convolution theorem (which is assumed throughout this paper) 


go 


dz t> 0. (8) 


Choosing T = 1 we have 


=m : 
t 


For any f(t) we have 


= f + lim2*, (10) 


t>0. (9) 


dp + tim” 


We thus have two ways of evaluating the z transform of f(#)/t.4 Wecan 
establish the above result quite independently by considering the con- 
volution integral in the p plane. In2z/(z — p) has a branch cut in the 
p plane extending from p = z to p = «. Since we are free to take the 
branch cut anywhere in the p plane, we will take it along the positive 
real axis. This is illustrated in Fig. 5. As before, the cross-hatched 


4 This equivalence should be utilized advantageously in evaluating the s transform of 
f(t)/t and also in summing infinite series. 
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area shows annulus of convergence and we have 


at 


wif 
x 
/ 


POLES OF 
p'F*(p) 


bRANCH CUT FROM 
p=: Z TO p = @ 


Fic. 5. Branch cut in p plane and path of integration. 


Along AB, p = |ple* ¢—0 
Along DE, p = | 


2 


Expansion for In — is given by (4) 
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Therefore 


2? 
= | (- 2) 13 
n n += + + (13) 


Now = 0. Hence, 
BCD 


J. (-2 


Now In k (where k is a complex number) = In |k| + arg k. It is 
then easy to see that Eq. 14 reduces to 


d 
p= 


One can generalize the above result for higher powers of ¢ in the de- 
nominator. We have 


p> nN 


Multiply by z—' and integrate from © to z. 


pr wap 


f(t) 
ie (16) 


Continuing this we obtain (1, 2) 


k times 


For k > 1 the above equation is defined for ¢ > 0; for k = 1,; Eq. 10 
is valid. 
f(t | 


form, actually it is not so because of difficulty in evaluating integrals 


Although the above results appear to yield z E fork > 1 inaclosed 


: 
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in a closed form of the type / (In (a + bx)/x). Therefore it is not 
possible to obtain closed expressions for z transforms of functions like 
f(t)/t* for k > 1.5 However, the usefulness of the result lies in evaluat- 
ing infinite series since numerical evaluations of definite integrals of the 


above type are available. 


3. zs TRANSFORMS OF FUNCTIONS SATISFYING 


Theorem: lf 


f(t, a) 
t 


eas = 


f F*(z, B)dB = ; f adds + lim? (19) 


n 


Note that F,*(z, 8), F*(z, a) are just z transforms of f;(¢, 8) and f(t, a) 
with a, 8 as parameters and are not to be confused with modified z 


transforms (2). 


Assumptions: 
1. f:(t, a) = 0 for t <0 which assures that F,*(z, 8) can be ex- 
panded in powers of z~' only. 


2. f(t, a) = 0 ast — 0 and lim— exists. 
no 
The importance of the assumptions becomes clear during the course of 


the proof. 
Proof: 


Taking the z transform on both sides (6) 
i, 
(20) 
n=0L amo ot 


Interchanging summation and integration with respect to a in Eq. 20 
(provided the integral exists) 


=0 


5 For k negative integer, the z transform is easily obtainable as explained in (2, 5). 


25 
: 
f(t, a) 18 
B)dp = ——. (18) 
t : 
: 
then 
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The term inside brackets on the left hand side of Eq. 21 is the z trans- 
form of B) 


a) 4 4, a) 
Fite, = tim 


nT 


We recognize that 


p-tDdp. 


The right hand side of Eq. 22 reduces to 


mi a) f(T, a)pdp +f f(2T, a)p-*dp 


(24) 


which is 


n—0 


f(0, a) = 0 by assumption. Hence, 


a 


Example: 


Consider the z transform of (1 — cos at)/t. Now 


0 1 0 t t 


Therefore 


1 — cos at 


pF*(p, a)dp + lim” (28) 


(1 — cos a)z (g + 1) 
(s — 1)(2? — 2s cosa + 1)’ 


(29) 


F*(z,a) = (1 — cosat) = 


(23) 
n 
: 
. 
where 
4 
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Instead of evaluating the integral with respect to z, we can evaluate 


f F,*(z, B)dB, where 
(z, B) 


zsin 
22 — 2zcos8 +1 


Since lim — 


z sin B be (1 — cosa)p(p + 1) 
[% (z? — 2z cos B + 1) ap - z (p _ 1) (p? Sie 2p cosa + 1) dp. (31) 


By actual computation these can be proved to be equal and therefore 


1 — cos at 2? — 2zcosa + 1 

————- = 3 In 32 

t 2—22+1 (32) 

It is seen that, in this case, integration with respect to a is slightly 
easier than integration with respect to z. 


Now 


s-* = $1n 


a=! n 2—22+1 (33) 


By manipulating this expression and substituting appropriate values of 
z and awe can sum upa number of infinite series. First puttingz = — 1 
we have 


(— 1)" = Incos “3 (33a) 
n=1 n 2 


Now integration with respect to a is definitely quicker to evaluate in- 
finite series since, as we see below, substitution of z = (— 1) in the 
integrand of the left hand side of Eq. 31 gives the summation very 
quickly. 
sin B — 1 
tan — dB = In cos B/2| 
= tn cosa/2. (336) 


This avoids the long process in evaluating the definite integral 


(1 — cos a)p(p + 1) 
@ 


Let us derive a few well-known series using the result 


1 — cosan — 2zcosa + 1 
n 2 —22+1 


; 
j 
; 
(30) 
; 
ie 
. 
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If z = 1 then the series diverges on the right hand side. 
Putting z = (— 1) in Eq. 34, 


1)" = In COS 5. 


When a = 2/2, 


that is, 


1 1 


When a = 


that is, 
—2(1+3 . (36) 


Therefore the sum ™ (1/(2n + 1)) does not converge. Hence, by 
n=1 


giving any value of a we can get a closed form and it also tells us whether 
the series converges or diverges. Thus Eq. 26 furnishes us two ways 
to evaluate an integral providing (18) is satisfied and as the example 
above has shown, we can use the equivalence of the integrals to sum 
up certain infinite series very easily. 


4. SUMMATION OF INFINITE SERIES 


One of the problems in the study of infinite series is to obtain closed 
expressions and operational methods have been used to obtain such 
results (7, 8). 2 transform by its very definition lends itself very easily 
to such summations and there are extensive tables of z transforms for 
various discrete functions of time. In this section we will use the re- 
sults of the previous sections to develop closed expressions for summing 


certain types of infinite series of the type } (f,/n*).6 Summation of 
n=O 


* The summation of the type Py (f(n + m)/n*),0 < m < 1, can also be obtained by sim- 
ple modification of Eq. 17 to read: 


k times 


F*(p, m) = 23mf(t)| 
3nf(t) = modified (2) z transform of f(t). 


n=] 
1 — cos=n 

| —cosnr 

—— 1)" = Incos= = ~, 

where 
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series of the type ¥ fag, using convolution z transform has been dis- 
n=0 


cussed in the literature (9, 10). 


Example 1 


Using convolution z transform we have developed 


0. 
t> 


Substituting z 


—1 


(— = In2. (38) 


n=l 


Example 2 


Letting 2-' = y in Eq. 37 and integrating with respect to y from 
JS, we obtain 


vo] y 1 
48 In - 


1 
=(1-—y)n(i-y)+y 


= (1 — 27) In (1 — + 2. 


With z 
1 
+ 1) is 
With z (— 1), 


Example 3 


Multiplying by z-' and integrating between the limits /,° we get 


nai 7? p —1 


3 
and 
1. (40) 
1)*+ 2in 2 — }. (41) 
a n(n + 1) : 
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This integral is not immediately evaluable. However, changing the 
variable of integration on the right-hand side from p-! = x reduces the 
integral to the form fo?‘ x-' In (1/(1 — x))dx. Tables are available for 
this integral between the limits /o'. Substituting z = 1 and using 


tables, we have 
Tr 


(43) 


6° 


Summation > (1/n*) for even values of k > 2 has been developed using 


other techniques (4, 7, 8). 


Example 4 


Evaluations of the following infinite series can be done easily by 
manipulation of Eq. 37. 


= (n+ I)(n+2) _ (44) 


(45) 


(46) 


Example 5 


To evaluate >> | ——— 
n 


n=0 


3 (47) 


n=0 n 


Recognizing that 


= cos BtdB, (48) 


we can evaluate the z transform in either of the two ways mentioned in 
Sec. 3. If we do that we get 
(49) 


a+-—lIn 


ete 
For a = x/2 andz = 1, 


2 


n=0 nN 


: 

| 

n=0 

> n+ 1 4 = 

| 

= : - 

if 

3 

5 + — In e?(*/2) 

3a 

= —, (50) 

4 

: 
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Recognizing that 


(52) 


This series known after Leibniz (11) is known for its clarity and sim- 
plicity and also affords a way of determining the value of 7. 


Example 6 
The z transform for (f(t)/(¢ + a)) has been derived (2): 


where a is a positive constant > 0. 


If f(t) is a unit step function then F*(z) = = i and assuming a to 


be positive integer 


This result can be easily verified by substituting the known series for 


In a and getting back the result > (1/(n + a))z-". 
n=0 


; 
sin 
2 us 
ae lim ——— = =, (51) 
nN 2 
we get 
1 1 
= 
: 
1 1 1 d 
+ + p 
po p p — | 
& 
—(a—1) —(a—2) —1 
[in 1 {2 1 a 2 + 1 (55) 
a—1 ] : 
| p= P 
5 
- 
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Example 7 
z transforms of (1/(m + a)(n + 6)), with a, b positive integers, are 


1 o g-n 1 1 
(56) 


As an example let a = 1, b = 2: 


1 2 


— g2 
z+ (zg 2*) In 


+ 1-2) |. (58) 
Similarly we can develop 4 transforms for (1/(m+a)(n+b)(n+c)---). 
The above results are in agreement with those of Stone (12). 


If we substitute z = 1 we have from (58) 


x 1 


With z = — 1 we have 


(— 1)" 
(n + 1)(n + 2) 


= 2In2 — 1. (60) 


The above examples have shown that summation of infinite series is 
equivalent to evaluating a certain integral and by obtaining closed 


forms for expressions of the type > (f,/n*)z-" we are able to tell 
n=O 


whether a particular series is convergent or divergent. 


CONCLUSIONS 


In this paper have been shown: (a) the equivalence of two forms of 
integrals, which could be used in evaluating an infinite series in a simple 
manner; (6) the problems associated with and care required in inte- 
grating functions involving branch point in the z plane. The integra- 


rt 
g 
i 
= 
‘ 
: 
: 
4 
+ 
: 
i 
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tion is useful in obtaining the z transform of functions with an irrational 
power of z. Further investigation into evaluating such functions is 
indeed warranted. Finally, (c) a closed form of infinite series using the 
convolution integral is obtained and a few examples are discussed. 
Extension of this method to obtain closed forms for certain other series 
is possible and conditions on convergence or the possibility of obtaining 
a closed form of the series is readily assessed from the form of the inte- 
grand in the convolution integral. Further extension of the convolu- 
tion for obtaining a closed form of a finite series of the form 


is possible and worthy of further investigation. 
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DESIGN RELATIONS FOR RESONANT POST 
WAVEGUIDE FILTERS 


BY 
ROBERT D. WANSELOW ' 


SUMMARY 


This paper presents a simple design procedure for resonant post waveguide filters. 
Such microwave circuits commonly are physically realized as ladder networks con- 
sisting of cascaded single post resonant cavities or double post resonant elements 
coupled by quarter-wave sections of transmission line. With this simplified tech- 
nique, design effort is minimized by the use of pertinent parameter relationships 
expressed in the form of a few graphs. A comparison is made of the theoretical and 
measured characteristics of two experimental filters each employing one of the two 
post configurations. The results show that this design procedure is quite adequate 
for most engineering applications. 


INTRODUCTION 


This paper considers waveguide filters that are of the band-pass type. 
These filters employ inductive posts as the primary reactive element in 
either of two basic configurations. The first and most common con- 
figuration is the cavity type where two identical single posts (longi- 
tudinally spaced approximately a half wavelength along the guide) 
make up one stage of the filter. The other configuration is the resonant 
element type where one filter section consists of a pair of identical posts 
(double posts) positioned transversely in the guide. This configuration 
is less known, but exhibits many desirable characteristics. Multi- 
stage filters are considered where both of these primary post configura- 
tions are employed as cascaded resonators coupled by quarter wave- 
lengths of transmission line as shown in Figs. la and 1b. The theory of 
quarter-wavelength coupling has been discussed by many (1, 2, 3)? quite 
adequately. As the Q of a given microwave filter is increased, the 
fabrication and adjustment of directly coupled resonant structures 
become increasingly more difficult. Therefore, by employing quarter- 
wave coupling between the filter resonators, the mechanical tolerances 
on the dimensions of the coupled elements are eased with no noticeable 
degradation of filter performance. This is especially true where single 
post cavity-type filters are required with loaded Q’s ranging:from 15 to 
200. Since the double post resonant structure is inherently a low Q 
type circuit, the quarter-wave coupling of this type of resonator pro- 
duces a filter which is reduced in length by a factor of two over a direct 
coupled cavity filter and by a factor of approximately three over a 


1 Microwave Tube Division, Hughes Aircraft Company, Los Angeles, Calif. 
? The boldface numbers in parentheses refer to the references appended to this paper. 
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quarter-wave coupled cavity-type filter, each possessing the same 
number of sections. The double post configuration is a very applicable 
type resonant circuit for filters whose loaded Q’s vary between 4 and 20. 
It should be noted that the correlation between theory and practice 
may be in error by as much as 5 per cent when either of these basic 
post configurations is employed when requiring loaded Q’s near their 
physically realizable limits. 


Pai Tuning Screws 


a, n - Section Quarter - Wave Coupled 
Single Post Resonant Cavity Filter 


Screws 


ty 
hd 


b, n- Section Quarter- Wave Coupled 
Double Post Resonant Element Filter 


Fic. 1. Two principal types of quarter-wave coupled resonant post waveguide filters. 


The tuning of a single post resonant cavity is accomplished by one 
tuning screw symmetrically placed between the two identical posts of 
the cavity, and with a double post resonant element a screw is also 
placed between these two transversely positioned posts. Mumford (4) 
and Reed (5) have discussed to some extent filters employing single post 
resonant cavities and Alstadter and Houseman (6) have suggested the 
use of double post filters. The theory of the double post as a wave- 
guide obstacle was developed by Gruenberg (7). Thus this paper 
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coordinates some of this information previously published together 
with some additional experimental design data which have proven very 
valuable for the rapid design of waveguide filters employing posts as 
the basic reactive waveguide obstacle. 

The simple design procedure presented here details the development 
of quarter-wave coupled waveguide filters employing either of the two 
basic post configurations described. Very useful experimental curves 
of the loaded Q, Q1, vs. physical dimension ratios of the posts are given. 
To the author’s knowledge, this information has not previously been 
published for the double post resonant element configuration. 


Fic. 2. General forms of low pass prototype filters. 


DESIGN PROCEDURE 


The first step in the synthesis of band-pass filters is the selection of 
the transmission coefficient corresponding to a network having the 
desired rejection-band characteristics. Two of the commonly used 
transfer functions are the Butterworth or maximally flat and the 
Tschebyscheff distributions. For given band-rejection characteristics 
or skirt selectivity, fewer resonant elements or stages are required when 
the synthesis is based on a Tschebyscheff distribution of the first kind 
with small tolerable ripple rather than a Butterworth distribution. 
Thus, upon selection of an appropriate transmission coefficient, the 
corresponding low pass prototype -section filter, as shown in Fig. 2, 
can be synthesized from which a low-pass to band-pass transformation 
is to be made. From the low-pass prototype the normalized reactive 
parameters, g;, can be obtained for use in the following design procedure. 
The transformation from a low-pass ladder network to its transmission- 
line quarter-wave coupled equivalent circuit of a band-pass filter is 


1 g 2 g 
i 
| | 3 3n-1 
at 
n even 
1 92 Gn-] 
g 93 Gn 
«= 
n odd 
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shown schematically in Fig. 3. The loaded Q, Q:x, of each resonator 
will be defined in the formulation of the following design equations. 
The effective Q, Q:, of the k* resonator of either a quarter-wave 
coupled cavity or resonant element type band-pass filter is equal to the 
loaded Q, Q:i, of the k* resonator, plus the Q due to the quarter-wave 


° Ago 

2 ge 


| 


| 
| 


Fic. 3. Quarter-wave coupled equivalent circuit of a band-pass filter. 


Ly. 


lines, Q,, which are connected to this resonator. Thus, we have the 
following relation : 


Qe = Que + Qa (1) 
and, by transposing, 

Que = Qk — Qa. 
The loaded Q, Qo, of the filter may be stated as follows: 


(3) 


where wo is the center frequency and w,’ is the 3-db bandwidth of the 
band-pass filter. The effective Q, Q:, of the kt" resonator may be 
equated as follows: 

(4) 
where g; is the normalized reactive low pass prototype parameter and 
w’ is the ripple tolerance bandwidth for Tschebyscheff designs and the 
3-db bandwidth for maximally flat or Butterworth distributions. 
From Cohn (8) a relation between normalized values of w,;’ and w’ is 
shown in Eq. 5, 


3 db = 10 logio [1 + (104° — 1) cosh? cosh a) (5) 
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where A,, is the ripple tolerance level in db corresponding to the normal- 
ized w’ = 1 point, m is the number of resonators in the filter, and 
a= =. Some bandwidth correction curves for converting ripple 


tolerance band width to 3 db bandwidth are shown in Fig. 4 from which 


Tolerance 


Am, 


Fic. 4. Bandwidth correction curves. 


a (the bandwidth correction factor) may be obtained. An approximate 
form of the loading upon a given resonator by its respective quarter- 


wave connecting line is 
Xx 2 
6) 


where m is the integer 1 or 2 depending upon whether the resonator is 
connected to one or two quarter-wave lines, respectively, and go is the 
guide wavelength and \, the free space wavelength corresponding to the 
filter center frequency. Therefore, by combining Eqs. 3, 4 and 6, 
Eq. 2 may be summarized in a more useful form with equal input and 
output impedances 


(7) 


1 fork = landn 
2 fork = 2,3,...,n — 1. 


m = 
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Design curves for the loaded Q, Q1, vs. the d/a ratio for 
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inductive susceptances. 


4-4-4 


The relation between the loaded Q, Q1, and the normali 


susceptance, B, of a resonant cavity formed by two shunt} 
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It should be noted that a = 1 for a Butterworth distribution. Values 
of gs may be determined by synthesizing the desired transmission coeff- 
cient or by utilizing Weinberg’s tables (9). 

As developed by Reed, the relationship between the loaded Q, Q1, 
and the normalized shunt inductive susceptance, B, of a post employed 
as a cavity element is 


=1(™)'|- (tan 5) + (8) 


and is plotted in Fig. 5. Once the susceptance of a single post is known 
as a function of its diameter (10), then the design curves of Fig. 6 can 
be employed quite effectively. In all these graphs the ratio of f/f, is 
the ratio of the operating center frequency of the filter to the waveguide 
cut-off frequency. 

The length dimensions between posts of an n-section quarter-wave 
coupled single post resonant cavity filter are shown in Fig. la. These 
dimensions are readily computed from the following : 


Ly = — tan | (9) 


1,2, ..., m for the cavities and 


where k = 2, 3, ..., m for the quarter-wavelength coupling lines. The 
susceptance term, B;, can be obtained by solving Eq. 8 for the loaded Q 
of the k** resonator which has been graphically done in Fig. 5. The 
guide wavelength \,o at the center frequency of the filter is readily 
obtained from the following : 


Ago 3 (Agi + (11) 


where \,; is the guide wavelength at the lower 3-db band edge and )y4e 
corresponds to the upper 3-db frequency point. 

For quarter-wave coupled resonant element filters (double posts) 
experimental design curves of the loaded Q as a function of the posts’ 
transverse position is given in Fig. 7. Once the values for the loaded Q 
of each resonant section are obtained again from Eq. 7, it is only 
necessary to select the post dimensions from this family of curves and 
to space the pairs of posts by \,0/4 lengths of line as shown in Fig. 10. 
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Example 


Suppose one is required to design a three-stage quarter-wave coupled 
single post resonant cavity-type filter in X-band waveguide with the 
following basic specifications: 


Center frequency = 9100 mc 

3-db bandwidth = 100 mc 

Ripple tolerance 1 db is permissible 
Number of stages = 3 


Procedure 


1. From Weinberg’s tables 
prototype parameters are: 


(9) the normalized reactive low-pass 


2.0236 
0.9941 
= g, = 2.0236 


3. From the bandwidth correction curves in Fig. 4, a 3-section filter 
(n = 3) with a 1-db ripple tolerance level (A,, = 1 db) exhibits a value 
a of 1.094. 

4. The center guide wavelength, d,o, can be obtained from Eq. 11 


where 


Ay: = 1.892 in. for f; = 9050 mc 
Age = 1.849 in. for f2 = 9150 mc 


in standard RG-52/U guide. Thus: 


Ago = 3(1.892 + 1.849) 


1.870 in. 


which corresponds to a free space wavelength, Xo, of 1.296 in. 


103 
ees 2. The loaded f the filter is, f Eq. 3 Aa 
e loade , Yo, of the filter is, from Eq. 9, 
2 9100 
Qo = = 91 ‘ 
=~ x = 
2x 100 
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5. From Eq. 
calculated. 


7 the loaded Q, Qi, of each resonator can now be 


2.0236 1.870 \? x 

Qu = 1.004 ( 2 ) 8 
= 99,7 


47.8 


Ors = 99.7 


6. The ratio of the center frequency to the waveguide cut-off fre- 
quency is 


i, 


=) 


7. Thus from the curves of Fig. 5, using f = 1.4, the normalized in- 


ductive susceptance can be obtained 


B, 
B, 


B; 


8. From the curves of Fig. 6 the : ratio of each post can be obtained 


0.132 


0.098 


= 0.132 
a 


a 


Therefore, since a = 0.900 in. 


d» 


d; = 0.119 in. 
0.088 in. 


4 
é 
: 
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7.8 
m= B, = 7.8 
a 
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= 
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9. From information obtained in steps 4 and 7, Eqs. 9 and 10 can 
now be employed to determine the space separation of the posts: 


= 0.860 in. 


1.870 
2X 
0.830 in. 
0.860 + 0.830 1.870 
2 4 
0.378 in. 


7.8 


— 
5.3 


The dimensioning of an n-section filter of this type is shown in Fig. 1a. 


EXPERIMENTAL RESULTS 


Figure 8 shows the measured and theoretical insertion loss response 
characteristics of a 3-section quarter-wave coupled single post resonant 


Insertion Loss in db 


Normalized w 


Fic. 8. Tschebyscheff response characteristics of a three section single post resonant cavity filter. 


cavity-type waveguide filter where the pass-band ripple tolerance was 
designed to be 1 db. These curves exhibit a close correlation between 
design and theory ; the step-by-step design for this filter is given in the 
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preceding section of this paper. A comparison between the experi- 
mental and theoretical results of a 4-section quarter-wave coupled 
double post resonant element type filter is shown in Fig. 9. This rela- 


Insertion Loss in db 


-1.6 -1.2 -0.8 -04 ° +04 +08 +1.6 
Normalized w 


Fic. 9. Butterworth response characteristics of a four section double post resonant element filter. 


‘tively low Q filter (Qo = 9.1) was designed to exhibit a Butterworth 
distribution. Both filters were made in RG-52/U waveguide and tuning 
was accomplished with 4—40 screws. 

As the Q of a given microwave filter is increased, the fabrication and 
adjustment of the resonant structures become increasingly more difficult 
if they are directly coupled. Thus, by employing quarter-wavelength 
coupling of posts rather than irises, the fabrication process is made 
much simpler and cheaper. For permanently tuned filters the com- 
plete integral assembly has been manufactured in aluminum with 
aluminum tuning screws. Upon tuning the filter and locking the tuning 
screws in place the filter is then dip-brazed, after which any part of the 
tuning screw mechanism protruding above the outside waveguide sur- 
fave is machined off, resulting in a smooth outside wall surface. 


TUNING PROCEDURE 


The tuning procedure for quarter-wave coupled filters involves a 
method of quarter-wave shifting the VSWR nulls on the input trans- 
mission line to the filter. A standing wave detector is inserted between 
a matched generator and the input side of the filter with a matched 
detector mount on the terminating side of the filter. Initially, all the 


: 
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resonant elements are detuned by inserting the tuning screws across 
the guide with the generator adjusted to the desired resonant center 
frequency. The positions P and P’ of two minima of the voltage 
standing wave are determined. The first cavity or resonant element 
is then tuned with the first tuning screw until one of the minima coin- 
cides with a point R, halfway between P and P’, that is, until the 
standing wave minima are shifted by \,0/4. The second resonant 
circuit is now tuned with the second screw until the minima shift back 
to their original points, P and P’. Thus, successive resonators are 
tuned by alternately adjusting the position of a standing wave minimum 
to occur at points P and R. The last resonator is adjusted so that 
maximum power will flow to the terminating detector. Upon comple- 
tion of this tuning procedure, the filter should be checked at the edges 
of the pass band to make sure that the over-all behavior of the filter 
is satisfactory. 


CONCLUSION 
Use of the procedures and graphs presented enables the microwave 


engineer easily to design waveguide filters with loaded Q’s ranging 
between 4 and 200. This procedure is sufficiently adequate for most 


engineering applications. Graphs of the loaded Q’s as a function of the 
physical dimension ratios of the posts have been given in a form that 
can be readily applied to the design of a waveguide filter in any standard 
full height waveguide where the a/b ratio is approximately two to one. 


REFERENCES 


(1) W. L. PrircHarp, “Quarter Wave Coupled Waveguide Filters,’ J. Appl. Phys., Vol. 18, 
pp. 862-872 (1947). 
(2) R. M. Fano anp A. W. Lawson, ‘Microwave Filters Using Quarter-Wave Couplings,’’ 
Proc. I.R.E., Vol. 35, pp. 1318-1323 (1947). 
(3) G. L. RaGan, ‘Microwave Transmission Circuits,’’ M.I.T. Radiation Laboratory Series, 
Vol. 9, New York, McGraw-Hill, 1948, pp. 677-706. 
(4) W. W. Mumrorp, “Maximally-Flat Filters in Waveguide,’ Bell System Tech. J., Vol. 27, 
pp. 684-713 (1948). 
(5) J. REED, “‘Low-Q Microwave Filters,” Proc. I.R.E., Vol. 38, pp. 793-796 (1950). 
(6) D. ALSTADTER AND E. O. HovusEMAn, ‘‘Some Notes on Strip Transmission Line and 
Waveguide Multiplexers,” JRE WESCON Convention Record, Part 1, pp. 54-69 (1958). 
(7) H. GruENBERG, “Symmetrically Placed Inductive Posts in Rectangular Waveguide,” 
Canadian J. Phys., Vol. 30, pp. 211-217 (1952). 
(8) S. B. Coun, ‘‘Direct-Coupled-Resonator Filters," Proc. I.R.E., Vol. 45, p. 188 (1957). 
(9) L. WetnBERG, ‘‘Modern Synthesis Network Design from Tables I-IV,” Electronic Design, 
Vol. 4, Nos. 18-21, Sept. 15, Oct. 1, Oct. 15, November 1, 1956. 
(10) N. Marcuvitz, ‘‘Waveguide Handbook,”’ M.I.T. Radiation Laboratory Series, Vol. 10, 
New York, McGraw-Hill, 1951, p. 262. 


| 
| 
> 
| 

| 

| 

| 
it 


THEORY OF NONLINEAR CONTROL * 
BY 
¥. 


ABSTRACT 


In the introduction the author discusses the nonlinearity in Nature as exempli- 
fied by a simple pendulum, inverse square laws of gravitation and electromagnetism,’ 
electronics, and planetury motion. Nonlinear feedback control systems have been 
analyzed in symposia sponsored by Polytechnic Institute of Brooklyn, American 
Society of Mechanical Engineers, and Institute of Radio Engineers. Even nonlinear 
system synthesis has been attempted. 

Part II discusses the three stages in the study of physical systems: (1) Develop- 
ment of linear system analysis and idealization of actual systems as linear; (2) Recog- 
nition of nonlinear systems and analysis of such systems as they actually are; and (3) 
Introduction of the new concept of nonlinear control and improvement of the per- 
formance of linear and nonlinear systems. Under (3) it is recognized that besides 
linear compensation of systems with nonlinearities the logical alternative might be 
the nonlinear feedback control of physical systems for either stability or optimum 
response. 

Part III starts with a review of linearizing techniques developed by Kochenburger, 
Johnson, Klotter, and others, known as the describing function method. Other 
quasi-linearization, equivalent linearization, and piecewise linearization techniques 
by Chen, Booton, Bass, and Stern are discussed. 

Part IV deals with the phase-space, the switching surface, and related concepts. 
The method of simultaneous phase-plane equations in a multi-degree-of-freedom 
system is especially adaptable to nonlinear feedback control system analysis. The 
phase-space approach has been used by Kalman, Hopkin, and others for the design of 
nonlinear servomechanisms. 

Part V covers new analytic techniques and new transforms. Papers by Pipes, 
Madwed, L. J. Lewis, Gillies, Stout, E. Weber, and others are reviewed. Taylor- 
Cauchy transforms for nonlinear systems are presented. 

Part VI discusses theory of nonlinear systems, including recent developments in 
stability theory and criteria, optimum nonlinear control systems, sampled-data 
systems, and nonlinear systems with random inputs. These recent developments give 
substantial support to the author's statement: ‘‘The newest aspect of feedback control 
is the development of the theory of nonlinear control.” 


I. INTRODUCTION 
Nature Is Nonlinear 


As nature embodies a manifold of beautiful phenomena, it cannot be 
limited by a linear relationship between elementary quantities. Thusa 


* Presented at the First International Congress of IFAC on Automatic Control, Moscow, 
U.S.S.R., June 27 to July 6, 1960. The complete English version appears here for the first 
time. While a Russian translation of the complete paper is available at the Moscow Congress, 
an abridged English version is to appear in the IFAC Congress Proceedings, published by 
Butterworths Publications, Ltd., London. 

! The Moore School of Electrical Engineering, University of Pennsylvania, Philadelphia, Pa. 
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simple pendulum as discovered by Galileo represents a simple harmonic 
motion, and yet the equation of a pendulum 


6” + (g/l) sin = 0 (1) 


is a second-order nonlinear differential equation, in which @ denotes the 

angular displacement, 6’’ denotes d*@/di?, the angular acceleration, g 

denotes a constant due to gravity, / denotes length, and ¢ denotes time. 
Kepler’s laws of planetary motion can be expressed as follows; 


M 


Mie? — 0 (2) 


+ = 0 (3) 


where M, denotes the mass of a planet, 1 denotes the mass of Earth, y 
denotes the gravitational constant, r denotes the distance between the 
planet and Earth, and 6 denotes the angle. Both Eqs. 2 and 3 are 
nonlinear. It is of interest to note that not only Newton’s law of 
gravitation is an inverse square law, but also Coulomb’s law of force 
between electric charges and the similar law of force between magnetic 
poles are inverse square laws. 
Van der Pol’s equation 


’— —x*)x’+x=0 (4) 


representing an electronic oscillator, in which x denotes a normalized 
value of grid voltage (ratio of grid voltage to saturation voltage), 
is nonlinear and gives rise to limit cycles. Transistors and magnetic 
amplifiers are both nonlinear devices. 

It has been remarked that the tremendous progress in the physical 
sciences is in large part due to the fact that most natural phenomena 
are at least amendable to a linear approximation. So long as one care- 
fully notes the natural limitations of a linear spring (for which Hooke’s 
law holds) or the linear portion of a B—H curve (before saturation), 
simple relations can be utilized to study natural phenomena. However, 
one is at the same time aware of the complexity of nature and wishes to 
study nature as it really is. 


2. Feedback Control Systems 


It was only twenty-five years ago that H. L. Hazen published his 
paper entitled ‘“Theory of Servomechanisms” (1).2_ At Cranfield Con- 
ference 1951, G. S. Brown discussed ‘“‘Feedback System Engineering—A 
Challenging Educational Objective’ (2). Control systems are, in 
general, classified as open-loop (open-cycle) systems (Fig. 1) or as 


"2 The boldface numbers in parentheses refer to the references appended to this paper. 
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closed-loop (closed-cycle) systems (Fig. 2) depending upon whether 
a comparison of the controlled quantity with the reference or desired 
value is absent or present. Thus a system of time-operated traffic 
lights represents an open-loop regulator and a thermostat temperature 
control system represents a closed-loop regulator. According to Brown 
and Campbell (3), a servomechanism is defined as ‘‘an error-sensitive, 
follow-up, amplifying system permitting a wide range of input command 
remotely located from the element being controlled.” 

Truxal (4) pointed out: ‘The situations of primary interest to the 
designer of feedback systems are those in which he is using feedback for 
a definite purpose.’ This definite purpose falls into one of the five 


categories; 


1. Simplification of the practical realization of a specified transfer 
function. 
2. Reduction of the sensitivity of the system characteristics to 
changes in parameters. 
. Control over the impedance seen looking into the system at 


various points. 

. Simultaneous control over a number of different transmittances, 
corresponding to multiple inputs to the single system. 

. Control over the stability, or instability, of the system, as in an 
oscillator circuit. 


G(s) | 


Fic. 1. Block diagram of an Fic. 2. Block diagram of a 
open-loop system. closed-loop system. 


In the first application, feedback is used to overcome limitations of the 
physical components. In the other applications, feedback is intro- 
duced to effect a specific change in the characteristics of the system. 
The above remarks apply to both linear and nonlinear feedback control 
systems. 

N. Wiener (5) has borrowed the term cybernetics from Ampére to 
mean the science of ‘“‘control and communication in the animal and the 
machine.”” Thus a cybernetic system may refer to a system which is 
designed for the communication of information or to a system which is 
designed for the control of one or more variables. H. S. Tsien (6) 
remarked in his Preface to Engineering Cybernetics: ‘‘A distinguishing 
feature of this new science is the total absence of considerations of 
energy, heat, and efficiency, which are so important in other natural 
sciences. In fact, the primary concern of cybernetics is on the qualita- 
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tive aspects of the interrelations among the various components of a 
system and the synthetic behavior of the complete mechanism.” 
Wiener’s recent book (7) has examined the role of nonlinear process in 
physics, mathematics, electrical engineering, physiology, and com- 
munication theory. 

The emphasis on, nonlinear circuit analysis in automatic control 
systems in recent years can be seen from the following list of symposia : 


1. Conference on Automatic Control, Cranfield, 1951. 

2. Symposia on Nonlinear Circuit Analysis 1953 and 1956 (Poly- 
technic Institute of Brooklyn, New York). 

3. Symposium on Nonlinear Control Systems 1956 (American So- 
ciety of Mechanical Engineers). 

4. Fachtagung Regelungstechnik Heidelberg 1956. 

5. Symposium on Nonlinear Control 1958 (Institute of Radio 
Engineers). 


Il. THREE STAGES IN THE STUDY OF PHYSICAL SYSTEMS 


1. Development of Linear System Analysis and Idealization of Actual 
Systems as Linear 


The classical method of solving linear differential equations is well 
known and furnishes a great urge to the idealization of actual physical 
systems as linear. Thus in a mechanical system with a mass and a 
spring, a linear differential equation of the second order can be formu- 
lated, with M denoting mass, B denoting friction, and K denoting 
“spring constant.’”’ We have 


Mx" + Bx’ + Kx =f (5) 
where x denotes displacement, x’ = dx/dt denotes velocity, x’ = d*x/dt? 
denotes acceleration, and f denotes the forcing function or the applied 
force. Let the forcing function be a unit step, x(0) = 0, and x’(0) = 0. 


The solution is 


1 1 et 


where s; and Ss; are the roots of the characteristic equations Ms? + Bs 
+K=0. 

Heaviside’s operational calculus was developed by Carson (8), 
Bush (9), Berg (10), Cohen (11), Wagner (11a), and others (11b) in 
the years 1926-1942. From 1942 on, Laplace transforms (12) have 
been largely used in the study of feedback control systems. Thus, 
denoting s as a complex variable, X(s) as the transform of x(f), and 
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F(s) as the transform of f(t), Eq. 5 is transformed to 
(Ms? + Bs + K)X(s) = F(s) (7) 
giving 
F(s) + Bs + K (8) 


where G(s) denotes the transfer function of the system. 

Many techniques have been developed in the study of feedback 
control systems. To mention a few: Bode’s diagram (13), Nyquist 
criterion (14), Evans’ root-locus method (15), Mason’s signal flow 
graph (16). From 1947 on, control system synthesis rigorously calls for 
a logical procedure for the transition from specifications to system. 
The Guillemin-Truxal approach ‘(17, 18) proposed the following steps: 


1. The closed-loop transfer function is determined from the specifi- 
cations. 

2. The corresponding open-loop transfer function is found. 

3. The appropriate compensation networks are synthesized. 


It may be mentioned that Fourier transforms (19) and Fourier 
integrals (20) have been developed for the analysis and synthesis of 
networks. The frequency-analysis techniques are highly useful in the 
study of communication systems, such as modulation and random noise. 


2. Recognition of Nonlinear Systems and Analysis of Such Systems as 
They Actually Are 


Recognition of nonlinear systems as they are indeed has had a long 
and glorious history. Astronomers, physicists, and engineering scien- 
tists have attacked actual nonlinear problems in celestial mechanics, 
fluid dynamics, elasticity, electronic circuits, and many other fields 
with fruitful results. We shall roughly divide the methods of attack 


in four categories : 


1. Perturbation and iteration methods. 

2. Isocline method and other graphical methods. 

3. Linearization methods. 

4. Phase-space method, and the use of analog and digital computers. 


The perturbation method was developed by H. Poincaré (21) and 
A. Lindsteht (22) primarily for astronomical problems. Poincaré’s 
Les méthodes nouvelles de la mécanique céleste and Liapounov’s Probléme 
général de la stabilité du mouvement (23) are two classics which appeared 
in 1892-3. Lord Rayleigh’s Theory of Sound (24) was published in 
1893-96, in which nonlinear problems were discussed. 
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The term “‘isoclinal lines’ appeared in a paper by J. H. M. Wedder- 
burn (25) published in the Royal Society Proceedings at Edinburgh in 
1902. The isocline method was used by van der Pol (26) to solve 
Eq. 4 in 1921-6. In 1937 Andronow and Chaikin (27) published their 
Theory of Oscillations in which the phase plane and the method of 
isoclines were described. They referred to K. Runge’s Craphische 
Methoden (28) published in 1932 and remarked: ‘The best graphical 
procedure is the method of tsoclines,’’ although they were aware of 
Liénard’s construction (29) made known in 1928. The work of the 
Soviet Institute of Oscillations founded by L. I. Madelstam, of which 
Andronow and Chaikin were members, has become known to the 
world at large through N. Minorsky (30). Other graphical methods are 
developed by Ku (31), Hsia (32), Paynter (33), and Buland (34). 

The third category includes the work of Kryloff and Bogoliuboff 
(35) and the recent work of Goldfarb (36), Tustin (37), Oppelt (38), 
Kochenburger (39), Johnson (40) and others. Based upon Professor 
Lefschetz’s English translation, we know that Kryloff and Bogoliuboft 
considered primarily equations of the form 


x" + wx = ef (t, x, x’, €) (9) 


where «¢ is a small positive quantity and f is a power series in e, whose 
coefficients are polynomials in x, x’, sint, cost. (Asa matter of fact, 
generally f contains neither e nor ¢.) Quoting Lefschetz’s Introduction: 
“Similar equations are well known in astronomy and have been the 
object of systematic investigation by Lindstedt, Gyldén, Liapounoff 
and, above all, by Poincaré. In a general sense, one may say that the 
same methods are applied by Kryloff and Bogoliuboff. . . . The solu- 
tions are approximated by the first » terms of certain asymptotic 
representations; the first two terms usually suffice and yield what the 
authors call the ‘refined first approximation’ which they discuss at 
length.” 

In Chapters V and VI of reference 35, the method of linearization* 
is described, which frequently enables one to by-pass the differential 
equation and proceed directly from the physical problem to the ap- 
proximate solutions. Rewriting Eq. 9 as 


mx" + kx + ef(x, x’) = 0 (10) 


where m and k are positive constants, and obtaining as the first approxi- 
mation a solution x = a cosy, where a and y satisfy certain conditions, 
we may say that the first approximation satisfies to within the order é& 
the linear equation z 
mx" + dx’ + kx = 0 (11) 
3 See also Chapter XII of ref. 30 on ‘method of equivalent linearization” and “principle 
of harmonic balance.” 
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where, for » = Vk/m, 


d —* cos ¢, — av sin ¢) sin ¢ dd (12) 
avn Jo 


k 


k+ cos ¢, — av sin ¢) cos ¢ (13) 
0 


Thus the nonlinear restoring force term in Eq. 10 is replaced by \x’ +x, 
where ki = k —k. It can be remarked also that 5 = \/2m is the 
dissipation decrement in the equivalent linear circuit and w = Vk/m 
the proper period of its oscillations, to within the order of &. 

Kochenburger’s Describing Function Method (39) of analyzing 
nonlinear systems is based on three assumptions: 


1. There is only one nonlinear element in the system. 

2. The output of the nonlinear element depends only on the present 
value and past history of the input. No time-varying character- 
istics are included in the nonlinear element. 

3. If the input of the nonlinear element is a sinusoidal signal, only 
the fundamental component of the output is considered. 


The approximation involved in the describing-function method is 
similar to that used by Kryloff and Bogoliuboff. Independent de- 
velopments of essentially the same methods were made by L. C. Golfarb 
(36) in U.S.S.R., by Tustin (37) in Great Britain, and by W. Oppelt (38) 
in Germany. 

According to Tustin (37): 

It is evident that if the system does in fact maintain a self-sustained oscillation, the 
periodic variation of A (the input) must have a waveform such that, when distorted 

by the nonlinear relationship B = N(A) and when further operated upon by the re- 
mainder of the loop transmission Y(p), the waveform A is again produced, with the 
same fundamental and the same harmonics. 


Though this procedure of neglecting the effect of the harmonics 
produced by a nonlinearity cannot ‘‘be justified on any general grounds,”’ 
it may be justified in particular cases in consideration of the actual type 
of motion involved when self-sustained oscillation occurs. 

We shall review the phase-space method briefly. According to 
Andronow and Chaikin (27): 


Generally speaking a dynamical system with m degrees of freedom depends upon a 
PF certain number n of positional coordinates qi, ---, gn and the state of the system at 
time ¢ is fixed by the values of the g; and their velocities g;’._ We may consider the 
qi, gi’ as coordinates of a space S of 2m dimensions called the phase space. To each 
state of the system there corresponds the point M with coordinates q;, gi’ in the space 
S and as ¢ varies the point M describes a curve called a path, which describes the 
history of the system. . . . The totality of all the paths may be described as the 
phase portrait of the system. It represents all the possible histories of the system; 
any one of them is determined by a single state M,. . . .” 
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Although the phase-plane method is well known in the solution of 
second-order nonlinear differential equations, the actual solution of 
high-order nonlinear differential equations depends on the extension 
from the phase-plane method to the phase-space method (41-45). It 
may be pointed out that the phase-space ‘‘portrait’’ is a representation 
of a dynamic system instead of a graphical method. Either digital or 
analog computers may be used in connection with the phase-space 
method. 


3. Introduction of the New Concept of Nonlinear Control and Improvement 
of the Performance of Physical Systems 


The new concept of nonlinear control can be stated in general terms 
as follows (45): 


A feedback system can be controlled either linearly or nonlinearly. If the system 
has inherent nonlinear characteristics, it is possible to introduce linear compensating 
networks to improve the system performance. If optimum system performance 
desired in a control system cannot be obtained by a combination of linear transfer 
functions, it is most desirable to introduce nonlinearities into the control system. 
If the system has inherent nonlinear characteristics, linear compensation and non- 
linear control can both be introduced. In short, the newest aspect of feedback control 
is development of the theory of nonlinear control. . . . 


While a physical system with a given input has to be analyzed to 
know its response or synthesized to yield the desired output, a general 
nonlinear system with random inputs offers even greater challenge for 
analysis and synthesis. We shall here describe the Wiener-Lee-Bose 
approach (7, 46-48) briefly as follows. 

According to Wiener (7), the most general probe for the investiga- 
tion of nonlinear systems is gaussian noise with a flat power-density 
spectrum, because there is a finite probability that this noise will, at 
some time, approximate any given time function arbitrarily closely 
over any finite time interval. Now one can characterize a nonlinear 
system by its response to shot-noise excitation just as one characterizes 
a linear system by its response to a unit impulse or a unit step or a 
trigonometric input. The shot effect is to be characterized by a par- 
ticular set of orthogonal functions, namely, the Laguerre functions (46). 
In Fig. 3, the outputs of the Laguerre networks are denoted by the 
Laguerre coefficients “1, U2, The outputs have the follow- 
ing interesting properties: (1) they are gaussianly distributed; (2) they 
are statistically independent; and (3) they all have the same variance. 
Since the output y(t) of a nonlinear system depends on the past of the 
input x(t), which is characterized by the Laguerre coefficients, y(t) 
can be expressed as 


y(t) Flu, Ur, * (14) 
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View 
HERMIT 

SHOT X(t) |LAGUERRE| U2 Via) 

NOISE POLYNOMIAL Det 
GENERATOR NETWORK GENERATOR 
t PRODUCT A 
DEVICE (27) 
Fic. 3. Block diagram of the circuitry for the characterization of nonlinear systems. 


SHOT HERMITE 
U Vv 
NOISE X(t) _ | LAGUERRE —>| POLYNOMIAL 
GENERATOR NETWORK GENERATOR | (7) 
71 
AVE RAGING F(t) V (4) = 
| SvSTEM , DEVICE 
Fic. 4. Block diagram of the circuitry for the determination of the 
optimum system coefficients. 
Vio) GAIN 
ry 
HERMITE 4 
GENERATOR 
Ay 
UT EXPONENTIAL 
FUNCTION MULTIPLIER 
GENERATOR 
Fic. 5. Block diagram of the circuitry for the synthesis of nonlinear systems. 
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The Hermite functions are chosen as the expansion for the function F 
because they form a complete orthogonal set over the interval — © to 
«©, and are particularly adapted to a gaussian distribution. Equation 
14 can then be given in the form 


y(t) = lim > A.V (a)e“”? (15) 


where = + + + --+-, V(a) represents the output 
of Hermite polynomial generator, and A, represents the characteristic 
coefficient of the nonlinear system, as shown in Fig. 3. In the practical 
case, we use a finite number of Laguerre coefficients and Hermite 
functions. Thus the limit s — © may be left out. Using the least 
weighted mean square error criterion, we get‘ 


y(t)V(a) = Aa (16) 


which provides the basis for the experimental determination of the 
characteristic coefficients A,. In getting Eq. 16, we have replaced the 
time average by the corresponding ensemble average according to the 
ergodic theorem. 

The experimental procedure is shown in Fig. 3, as given by Bose 
(47, 48). Since y(¢), in an ideal case, should approach the desired out- 
put, z(t), the experiment can be carried out by taking the averages of 
z(t)V(a) instead of y(t)V(a). Thus we can have the procedure as 
shown in Fig. 4. After the characteristic coefficients A, are determined 
from the experiment shown in Fig. 4, a synthesis procedure can be 
developed according to Eq. 15 with s finite, as shown in Fig. 5. 

The Wiener-Lee-Bose procedure has been developed for the analysis 
and synthesis of optimum nonlinear filters. However, the method is 
applicable to the study of nonlinear control systems with random inputs. 
It is therefore possible to analyze or synthesize nonlinear control systems 
for random inputs and desired outputs, once one has decided to use non- 
linear control for the improvement of the performance of physical 


systems. 


Ill. LINEARIZATION TECHNIQUES 


1. Development of the Describing Function Method 


We have already mentioned the contributions of Goldfarb, Tustin, 
Oppelt, and Kochenburger in the development of the describing func- 
tion method. E. C. Johnson (40) applied the sinusoidal analysis to 
study the free play or backlash in servomechanisms. While in the 
contactor servomechanism Kochenburger chooses a linear compensating 
network, in the backlash case Johnson selects a suitable nonlinear 
mechanical linkage to fit into the system. 


4 y(t) V(a) denotes the time average of y(T) V(a). 
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Use of describing function representation of nonlinear elements 
substitutes a quasi-linear system whose parameters change slowly for 
one containing rapidly varying relationships. The new system with 
slowly changing parameters can be thought of in terms of the poles and 
zeros of its response functions. As amplitude and frequency shift, the 
poles and zeros of the hypothetical system move about in the complex 
plane. Thus stability of the system can be studied similarly as in the 
linear case. One can tell whether the hypothetical system is stable or 
unstable, or has periodic oscillation. However, as both Johnson and 
Nichols (49) pointed out, the low-frequency intersecting point of two 
loci is a divergent equilibrium point, while the high-frequency intersect- 
ing point is a convergent equilibrium point. The low-frequency di- 
vergent limit cycle prediction must be disregarded if the describing 
function method is to give predictions that check with exact analysis 
and actual physical experiments. 

At the Heidelberg Conference 1956, J. C. West (50) discussed the 
use of frequency response analysis in nonlinear control systems. K. 
Klotter (51, 52) at the same time offered a paper on Hamilton Describing 
Function in contrast to Kochenburger’s Fourier Describing Function. 
There have been so many papers using the describing function technique 
in recent years that we can say this method has been generally accepted 
as a useful tool in the study of nonlinear systems. 


2. Other Linearizing Techniques 


We shall briefly mention three examples: (a) quasi-linearization; 
(b) equivalent linearization; and (c) piecewise linearization. 

The quasi-linearization technique developed by K. Chen (53) 
can be used for transient study of nonlinear feedback control systems. 
Booton’s quasi-linearization method (54) for analysis of nonlinear 
systems with random inputs will be discussed later. Based upon the 
principle of quasi-linearization (53) that ‘‘the quasi-linear representa- 
tion of a nonlinear element should be obtained by testing the element 
with its actual input,’’ a nonlinear element for transient study should 
be obtained by testing the nonlinearity with a transient type of signal 
which simulates the actual input. Assuming that the nonlinearity is 
piecewise linear, the operating point of the nonlinear element is expected 
to shift during the transient over several linear ranges. A quick 
evaluation of the actual transient response in the first linear range 
through linear servo techniques determines the test function for the 
nonlinearity in the first two linear ranges. The quasi-linearization 
representation thus obtained permits a quick evaluation of the actual 
transient response in the second linear range. This will in turn de- 
termine the test function for the nonlinearity in the first three linear 
ranges, and so forth. 
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R. W. Bass (55) proposed the equivalent linearization which can be 
explained in a houristic way as follows. Consider the difference- 
differential equation 


Ax!" (t + ta) + (t + te) + Bx’ (t + ta) + Cx(t + ta) 
+ x’()] =0 (17) 


where F denotes a nonlinear function of x(t) and its derivative x’ (é), 
A, B, C are constants, and ¢, denotes a time-lead. Note that the orig- 
inal equation might have a nonlinear function of x(¢ — t,) and x’(t — ta), 
where f, denotes a time-lag. The steps are: (1) insert the complex 
solution x(t) = ae*‘ in Eq. 17 to obtain the formal “instantaneous 
characteristic polynomial’’; (2) replace complex numbers in the argu- 
ments of all nonlinear terms by their real parts, and weight each of these 
nonlinear terms by a factor of 2; (3) divide the polynomial by ae*'; 
(4) multiply it by d(vt) and average over a period (from 0 to 27); (5) 
separate this averaged polynomial into real and imaginary parts and 
find the roots a, v of the resultant pair of real equations. Then one can 
study whether each solution (a;, v;) is stable or unstable as in the studies 
of E. C. Johnson (40, 56). . 
Piecewise linearization has a long history. However, T. E. Stern 
(57) has developed piecewise-linear network analysis and synthesis 
with the view of bridging the gap between the very general studies of 
nonlinear systems, exemplified by the work of Wiener (7), Zadeh (58), 
and others. A systematic formulation of the problem leads to a suit- 
able symbolism and an algebra of inequalities (59). This formulation 
helps to broaden the scope of diode network synthesis in particular and 
of piecewise-linear synthesis in general. It is of interest to note that 
Mason and Zimmermann (60) have used this piecewise-linear method in 
the teaching of electronic circuits. 


IV. THE PHASE-SPACE METHOD AND RELATED CONCEPTS 


1. The Phase-Space Method as a General Formulation of the Nonlinear 
Physical System 


As mentioned before, a dynamical system with » degrees of freedom 
depends on n positional coordinates g;, and the state of the system at 
any time depends upon the values g; and their velocities g,’. Thus a 
phase-space of 2n dimensions is associated with the behavior of the 
dynamical system. By choosing a system of coordinates we do not 
necessarily have to solve the problem graphically. However, the usual 
association of the “‘isocline method”’ with the solution of a phase-plane 
equation suggests a close relationship between the phase portrait and 
the analytic solution. 
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In discussing K. Chen’s paper (53), R. E. Kalman said: 


It should be strongly emphasized here that the phase-plane, or rather, phase space 
is not just a name for a clever method directed toward solving particular problems, 
but it is the most fundamental concept of the entire theory of ordinary differential 
equations (linear or nonlinear, time-dependent or time-invariant) at the present time. 
This is because the phase-space is a way of representing (sometimes only in an abstract 
sense) all possible solutions of a differential equation. . . . 


2. Method of Simultaneous Phase-Plane Equations 


In ‘‘Analysis of Nonlinear Coupled Circuits—II” (43), the author 
developed the method of simultaneous phase-plane equations (see 
Section 10-6, ref. 45, pp. 234-243). For a system of three degrees of 
freedom, the phase space is characterized by six coordinates—three 
positional coordinates g,; and three velocity coordinates q;’. In an 
electric network of three loops, the phase-space coordinates represent 
three charges and three currents. In a coupled circuit of two loops, 
the phase-space coordinates are reduced to two charges and two cur- 
rents. Though this method is developed later than the method of 
high-order phase-space equation, it is logical from the viewpoint of the 
phase-space representation of the state of a general dynamical system. 

As shown in references 61 and 62, the method of simultaneous phase- 
plane equations is highly powerful in the analysis of feedback control 
systems with nonlinearity either in the forward branch or in the feed- 
back branch or in both. Taking a two-loop system for illustration, the 
phase-space equations are 


dq,’ 
18 
te, (18) 


dq.’ 
(19) 


where g,”’ and q,’’ denote respectively the second derivatives of g: and 
g2 with respect to time, and they are in general functions of q, and q,’. 


Fic. 6. Block diagram of a nonlinear control system. 


In a nonlinear control system shown in Fig. 6 (see Section 12-9, ref. 45, 
pp. 303-305), there are two nonlinear controllers Gy; and Gye in the 
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forward branch and the feedback branch respectively. The two vari- 
ables are e and c, corresponding to the positional coordinates q; and q2. 
The other two coordinates in the phase-space are e’ and c’, corresponding 
to q;’ and q,’. The differential equations are 


+ gaile’,e) = m(t) = Te!” +c” (20) 
c’ + = b(t) =r—e (21) 


where r denotes input or reference signal, c denotes output, e denotes 
actuating signal, b denotes the feedback signal, m, denotes the control 
effort, 7 denotes the time constant, and g,; and g,2 represent two non- 
linear functions. The c’”’ term on the righthand side of Eq. 20 can 
be eliminated by differentiating Eq. 21 with respect to time and sub- 
stituting in Eq. 20. Thus we have 


+ guile’,e) + Te’ +e (22) 


The phase-plane equations corresponding to Eqs. 18 and 19 are 


de’ _ f(r, r, cl, c) e) — —¢@ (23) 
de e 4 


dc (24) 


Note that c”’ in Eq. 23 can be expressed in terms of (r — e — gy2) just 
as in Eq. 24. 

The simultaneous solution of these equations calls for numerical 
techniques or the help of computers. The general formulation is exact 
and complete. The nonlinear functions can be given either in analytic 
form or in graphical form. It is hoped that the programming of phase- 
space equations will help the study of nonlinear control systems. 


3. Phase-Space Approach to the Design of Saturating Servomechanisms 


We shall cite two concrete examples about the phase-space approach 
to the design of servomechanisms. R. E. Kalman (63, 64) made a 
phase-plane analysis of automatic control systems with nonlinear gain 
elements and then developed a method for the design of second and 
higher order saturating servomechanisms. Kalman’s method is based 
on linear transformations in the phase space, corresponding with the 
partial-fraction expansion of transfer functions to separate natural fre- 
quencies, and makes use of the root-locus method for the qualitative 
study of close-loop stability. The analysis leads to replacing a high- 
order system with a second-order system which closely approximates 
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the former. The method is applied to the study and design of saturat- 
ing servomechanisms. The principle of analysis is based on the 
following : 


The trajectories of a nonlinear differential equation may be approximated with arbi- 
trary accuracy by breaking up the phase space into a sufficiently large number of 
regions, such that within each region the trajectories obey a linear differential equation 
with constant coefficients. This is equivalent to approximating nonlinear charac- 
teristic curves by straight-line segments. . . . 


An earlier paper by A. M. Hopkin (65) used the phase-plane approach 
in the compensation of saturating servomechanisms. A new device is 
developed and called an anticipator. It controls the application and 
reversal of the manipulated variable. It senses the point where the 
actual system phase-plane trajectory crosses the particular trajectory 
making up the critical boundary, that is, one of the two which go into 
the phase-plane origin. It there reverses the manipulated variable to 
bring the system operating point along the critical boundary to zero 
velocity with zero error. In this study the optimum performance is 
assumed to be that behavior in which the system returns to rest with 
zero Or minimum error in the shortest time following a step input. 
The anticipator unit utilizes a biased-diode nonlinear element. For a 
constant-velocity input, a narrow-proportional-band amplifier is sug- 
gested as the nonlinear network. 

A recent paper by Hopkin and Iwama (66) reported about the 
design of a predictor-type air-frame controller by phase-space analysis. 
As the method of calculating the phase-space trajectories from analytical 
expressions for e, e’, and b (the feedback signal to the predictor con- 
troller) is exceedingly tedious, most of the trajectories are obtained by 
the technique of ‘‘reversing the time axis’’ with the use of an analog 
computer. For this study made on an analog computer, a workable 
solution is found in the form of a cathode-ray oscillograph-photocell 
combination with a mask which changes shape as b varies to simulate the 
effect of the varying shape of the projections of the critical boundary 
surfaces with 6. The particular predictor is set up to handle the case 
where the controlled system is of third order with one highly oscillatory 
mode. Thus a control system of the optimum relay type is obtained to 
control the linear pitch acceleration of an air frame and to compensate 
automatically for parameter changes resulting from changing flight 
conditions. 


4. Switching Criteria for High Order Servomechanisms 


In the field of contactor or on-off servos, phase-plane techniques are 
used by MacColl (67) and Weiss (68) to determine the response and 
improve the stability of second-order servos possessing such non- 
linearities as finite torque reversal time and coulomb damping. D. C. 
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McDonald (69), A. M. Hopkin (65), N. S. Neiswander and R. H. 
MacNeal (70), however, used the phase-plane technique to optimize 
performance. In McDonald’s method a torque-reversing criterion is 
developed in which the sign of the torque applied to the load is de- 
termined by a nonlinear function of error and error rate. I. Bogner 
and L. F. Kazda (71) extended McDonald’s method to obtain unique 
switching criteria for contactor servomechanisms of order greater 
than two. 

Bogner and Kazda’s analysis presupposes a discontinuous forcing 
function, fixed in magnitude and reversible in sign, supplied to a 
dynamic system the operation of which is describable over every in- 
terval by a linear differential equation. This type of control, while 
nonlinear in the general sense of the term, is sectionally or piecewise 
linear. The method makes use of the particular phase space in which 
servo error is plotted against its higher derivatives. It is of interest 
to note that this phase space having the variable and its higher deriva- 
tives as the coordinates* corresponds to the phase space used in the 
author’s method of high-order phase space equation (see Section 10—5 
of ref. 45, and refs. 42-44). Let the third-order differential equation be 
given by 
+ 6 (25) 


where 6 denotes either position or velocity input, positive or negative 
in sign. The three phase-space coordinates are: u; = @, U2 = e’, and 
u; =e’. In discussing the switching criteria, Bogner and Kazda 
made use of the principal coordinate phase-space, where the new coordi- 
nates are: 0; = U, — U3, Ve = Uy, + U3, and v3 = uz. In an isometric 
drawing of the v-space trajectory for a third-order contractor servo, a 
representative point moves from 0 to 1, where it intersects the first 
switching surface. The sign of the controlled variable is then reversed, 
permitting the representative point to move in the surface until it 
reaches the second switching point 2. It then follows the trajectory 
into the origin. 

N. J. Rose (72) used a more complicated switching rule, with the 
correcting force depending on both the error and the error rate. Thus 
it is possible to find the optimum switching rule which reduces the error 
and the error rate simultaneously to zero in a minimum time. The 
optimum switching criteria given by Rose apply to second and higher 
order systems. 

S. S. L. Chang (73) discussed the optimum switching criteria for 
higher order contactor servo with interrupted circuits and asked the 
question whether the assumption that the error and its derivatives up 


5 This phase space is a consequence of Cauchy’s Theorem of Existence for the solution 
of'a differential equation, as a differential equation of the n*® order is equivalent to a system 
of m equations of the first order. 
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to the (n — 1)* order should be continuous at the moment of switching 
is necessary or not. He investigated the optimum open-circuit switch- 
ing criteria for a high-order contractor servo by an extension of the 
phase-space method. Hung and Chang (74) extended the optimum 
switching theory to a more general form such that phase space analysis 
is applicable to a system which has non-vanishing power of s in the 
numerator of its transfer function. They concluded: ‘‘The phase-space 
method is a basic way of thinking in analyzing nonlinear system.” 


Vv. NEW ANALYTIC TECHNIQUES AND NEW TRANSFORMS 


1. Reversion Method and Iteration Method 


L. A. Pipes has applied operational methods (75) to the analysis of 
nonlinear systems. The reversion method (76) is based on the algebraic 
procedure used in reverting power series. Let the nonlinear differential 
equation be written in the form 


ay tay? + --- +asy°+--- =k (26) 


where the a; coefficients are, in general, operators and functions of the 
operator D = d/dt. Assume a solution of the form 


y = Aik + Ack? + Ask? + (27) 


The A, coefficients may be determined by substituting Eq. 27 into 
Eq. 26 and equating coefficients of equal powers of k. This method 
can be compared with the method of Picard. Whereas the reversion 
method is used to analyze a series circuit with a nonlinear inductor, 
Pipes introduced the minimum mean square error method in the 
analysis of a series circuit with a nonlinear capacitor (77). It may be 
noted that the author used the acceleration-plane method to analyze a 
series circuit with both nonlinear inductance and nonlinear capaci- 
tance (78). 

Lalesco’s nonlinear integral equation was also explored by Pipes (79). 
In the study of forced oscillations of nonlinear circuits, the method of 
undetermined coefficients was extended (80). For circuits with time- 
varying parameters, the Brillouin-Wentzel-Kramers method was sug- 
gested (81). An iterative method (82) for nonlinear vibrations was 
given by R. E. Roberson who earlier contributed a paper on the syn- 
thesis of a nonlinear dynamic vibration absorber (82a). 


2. Number Series Transformation Method 


The number series transformation method was introduced by Tustin 
in 1947 (83). A. Madwed (84) made further developments and gave 
the following advantages of this method: (1) the problems involved in 
analyzing and synthesizing systems are solved in the time domain; 
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(2) the method permits the analysis and synthesis of nonlinear systems ; 
(3) the actual mathematical operations involved are simple addition, 
subtraction, multiplication and division of real numbers; and (4) the 
method is easy because it is an operational method. It is shown (84a) 
that it is possible to take the Fourier and Laplace transforms of the 
number series transformations of functions. 


3. Harmonic Analysis and Periodic Solutions 


In order to justify equivalent linearization, Kryloff and Bogoliuboff 
(35, 30) observe that a nonlinear oscillatory process is generally charac- 
terized by a certain Fourier spectrum of the component frequencies 
resulting from the nonlinearity of the system. The principle of har- 
monic balance is assumed a priori, so that the equivalent parameters can 
be determined in the linearized system which gives the same oscillation 
as the fundamental harmonic of the quasi-linear system. L. J. Lewis 
(85) suggested the use of polynomial functions for carrying out the 
harmonic analysis instead of the more conventional use of trigonometric 
functions. In discussing this paper, Stout remarked: ‘‘In addition to 
the suggested applications to nonlinear differential equations, modula- 
tion, and distortion, it should be useful for the calculation of ‘describing 
functions’ employed in nonlinear feed-back system analysis.’’ It may - 
be noted that Lewis’ polynomial functions are closely related to the 
Tschebyscheff polynomials. 

Given a modified form of the van der Pol equation 


x” — + — px?)x’ + x = 2Bw, sin wit (28) 


where ¢ and uw are parameters which are assumed to be small, and uy is of 
the order «. A. W. Gillies (86) discussed the periodic solutions of the 
above equation which differs from Eq. 4 by the presence of a forcing 
function and a relatively large term pkx in the damping coefficient. 
The effect of this is that the solution has to be developed in powers of 
u, and the equations of the first approximation cannot be obtained with- 
out considering the second harmonic. Consequently methods such as 
those of Kryloff and Bogoliuboff or Cartwight’s difference equation (87) 
have to be carried a stage farther than is necessary for the van der Pol 
equation. In an earlier paper (88), Gillies used the power series in the 
study of nonlinear circuits in which a triode is associated with a linear 
network of any form. Part 1 deals with a stable circuit having zero or 
negative feedback ; Part 2 deals with the free oscillations of an unstable 
regenerative circuit; and Part 3 deals with the forced oscillations of a 
regenerative triode circuit. It is shown that the mathematical tech- 
niques applied in Parts 1 and 2 can be applied in Part 3, and the uniform 
procedure leads to a physical description of the phenomena in terms of 
modulation products. 
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4. Step-by-Step Method 


Many experts on numerical analysis have worked on the step-by- 
step method. The high speed electronic digital computer uses the 
step-by-step method. Whether one gets solutions of differential equa- 
tions in the time domain or one solves phase-space equations, a step-by- 
step procedure can always be resorted to. This procedure is also used 
in many graphical methods. 

We shall mention as an example a step-by-step method developed by 
T. M. Stout (89) for transient analysis of feedback systems with one 
nonlinear element. The method is somewhat similar to methods 
described by Tustin (83, 90), Madwed (84), Wagner (91), and Truxal 
(92). Unlike these methods, Stout’s method does not require resolution 
of the variables into step, pulse, or other component signals, or the sub- 
stitution of difference for differential operators; no attempt is made to 
use regression formulas or z-transform techniques. The weighting 
function required in the analysis is determined by inverse Laplace trans- 
formation of the appropriate transfer function, and the superposition 
integral is evaluated by standard numerical methods. 


Fic. 7. Block diagram of a system with one nonlinear element. 


Referring to Fig. 7, where y(t) is a nonlinear function of x(t) and 
can be expressed as y = N(x), we have 


= lt) = — fou (29) 


Note that the weighting function h(t) is the impulse response of the 
feedback block and is the inverse Laplace transform of H(s). An 
approximate solution is found as follows: Using numerical methods, 
the various system variables are computed at intervals At = T; these 
values are denoted by v, = v(m7T), h, = h(nT), and so forth. If the 
trapezoidal rule is used, the superposition integral is approximated by 


If we let 8, = Th,, Eq. 30 can be regarded as the product of the two 
sequences yo/2, yi, and Bo/2, Bi, Bn. The first value of 


® See also Naumov’s method (ref. 89b) and Dietz’s method (ref. 89c). 
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y(t), yo, is found from the initial conditions; the rest of the values are 
determined step-by-step. 


5. Complex Convolution Method 


The method is developed by E. Weber (93). The Laplace transform 
of the product g(t)h(t) is given by the complex convolution integral 


= G(s — = G(s)"H(s)_— (31) 


cm yo 


where the last term is a symbolic notation. Similarly, we have 


L[h(t)h(t)] = H(s)*H(s) = H®(s) 


(32) 


and 


efhr(t)] = H(s). (33) 


Taking for example the differential equation 


+ ayy’ + + ky) = f(t) 


(34) 


and applying the Laplace transform, we obtain 


+0) Y(s) +aok (s) = F(s) + (s+a1)yo+yo' = Fi(s) 


(35) 


where Y(s) and F(s) are the Laplace transforms of y(t) and f(¢), re- 
spectively, s is the complex variable, a; and a» are constants, yo and yo’ 
are the initial values of y and y’, respectively. Assume “rather boldly”’ 
that Y(s) be a mermorphic function and write it 


For n = 3, we get 


It has been assumed that there exists the Laplace transform of y(t). 
If Y(s) is mermorphic, then the Laplace transform of any integer 
power of the function is again mermorphic. Substituting Eqs. 36 and 
37 in Eq. 35 gives, for F(s) = 1/s and n = 3, 


aok A.A 


_ 1 + s*yo + s(aryo + 
s(s — s’)(s — 


(38) 


} 
: 
: 
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where s’ and s” are the roots of the characteristic equation s? + ais 
+a) = 0. With further simplification of the triple summation term, 
Eq. 30 can be reduced to a form where coefficient comparison is possible, 
provided the unknown s,’s are not identical to s’ or s”’. 


6. Taylor-Cauchy Transform Method 


The Taylor-Cauchy transform method (94) was developed by Ku, 
Wolf, and Dietz. The Taylor-Cauchy transform is defined by’ 


where \ is a complex variable which replaces the real variable ¢, and 
where F(\) replaces any complicated time function, like D?y = y” and 
N(y) = y*, by substituting W for y and X for t, like D?W(\) = W” 
and N(W) = W*. Taking Eq. 34 as an illustration, the transform 
equation is 


+ aD + ac)W()] + TLakW*(d)] = TLF(A)] (40) 


where D = d/dd and D? = d?/d)d*. The new transform is linear in the 
sense that it has additivity and communitivity properties. In general, 
the highest derivative term W‘*)(\) is assumed to be analytic for a class 
of systems. Then 


Applying the Taylor-Cauchy transform to both members of the above 
equation gives 


TLW™(d)] = Wau. (42) 


The inverse Taylor-Cauchy transform is then given by 


= = WH). (43) 


n=0 


The transform of nonlinear terms is discussed in refs. 94, 95, and 96. 
Three tables of Taylor-Cauchy transforms are given at the end of ref. 94. 
Given the simple nonlinear differential equation 


dy 
dt 


+ y? = u(t) = unit step. (44) 


7 C denotes a closed contour in the \-plane enclosing the singularities of F(A). 
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Taking the Taylor-Cauchy transform gives 
Wh + (45) 


in which 
n—2 
Ws Wn—p—2 
Ca-2,0” = > 


(46) 


Solving Eqs. 45 and 46 for w, recursively gives wo = 1, w; = 0, 


WwW, = —1, w, = 0, ws = 2/3, ws = 0, we = — 17/45, and so forth. 
Since w,’s are solved for the first derivative W’(X) in this example, 
we obtain 
W,A"t! 
W(A) = + (47) 


For yo = 0 and hence Wy = 0, Eq. 47 becomes 
2 


3 15 315> tanh (48) 


which gives y(t) = tanh ¢ as the solution of Eq. 44. 


VI. THEORY OF NONLINEAR SYSTEMS 
1. New Concepts and Theorems 


Some general theorems for nonlinear systems possessing resistance 
or reactance were given by Millar (97) and Cherry (98) in 1951. Energy 
relations in nonlinear inductors and capacitors were discussed by Manley 
and Rowe (99) in 1956. In a paper presented at the Cranfield Con- 
ference (100), Cherry and Millar gave nine new theorems and explained 
the following new concepts. 

The concept of content and co-content. Associated with any non- 
linear resistor characteristic, the content G = /vdi and co-content 
J = fidv are duals. The sum (G + J) always equals the instan- 
taneous power. 

The concept of energy and co-energy. The energy stored in an 
inductor is 7 = fid¢ and the energy stored in a capacitor is U = /vdq. 
The co-energy for an inductor is 7’ = /¢di and the co-energy for a 
capacitor is U’ = /qdv. 

We shall quote two Theorems as follows: 


Theorem 3. In any all-resistor multi-mesh circuit driven by current 
sources only, the total content G is unaltered if infinitesimal increments 
be given to any of the currents in the system, provided these increments 
are consistent with Kirchhoff’s current law. Similarly, in a voltage- 
driven circuit, the total co-content is stationary with respect to con- 
sistent voltage variations. 
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Theorem 8. In any circuit of voltage (current)-driven nonlinear 
inductor elements, the total energy T (co-energy 7’) has a stationary 
value with respect to node-voltage (mesh-current) variations, provided 
that these are consistent with Kirchhoff’s voltage (current) law. 

A dual to Theorem 8 gives Theorem 9, which is obtained by writing 
capacitor for inductor, node for mesh, voltage for current, charge for 
flux. In the case of any nonlinear system, possessing dissipative and 
energy-storing elements, the appropriate Lagrangian function L may be 
shown to be either L = 7’ — Uor L = U’ — T, a mixture of energies 
and co-energies. In terms of current variables Lagrange’s equations 
involve the first form of L and the content G. In terms of voltage 
variables the equations involve the second form of L and the co-content 
J. When the networks themselves contain nonlinear elements, nothing 
resembling Thévenin’s theorem applies. However, the use of the 
terminal characteristic of an active nonlinear resistive network often 
simplifies graphical work. The new concepts and theorems hold not 
only for electrical networks, but also for mechanical and electro- 
mechanical systems. 

According to Manley and Rowe (99), the nonlinear reactive ele- 
ments may be regarded as modulators in which the input signal fre- 
quency and a local oscillator frequency combine to produce new fre- 
quencies. These are harmonics and sum and difference combinations 
of the two applied frequencies, one or more of these resulting sidebands 
being taken as the useful signal output. The results are thus applicable 
to nonlinear reactor modulators and demodulators and to magnetic 
and dielectric amplifiers. 

Let the characteristic of a nonlinear capacitor be given by v = f(g), 
where f(q) is single-valued. If the nonlinear capacitor is connected in 
a circuit containing two generators of frequencies f; and fo, in general 
all of the frequencies fn. = mf: + nfo will be present in the circuit, 
where m and n take on all integral values, positive, negative, and zero. 
Assuming that no other frequencies are produced by instability, the 
charge g, and hence the current 7 can be written as a double Fourier 
series, using both positive and negative frequencies. The voltage v 
across the nonlinear capacitor must consist of the same frequency com- 
ponents as the charge g. Thus v may also be written as a double 
Fourier series. 

Denoting the coefficients of the double Fourier series by Qa,n, Im.n; 
and V,,.,, respectively, for g, 7, and v, the average power associated with 
(mf, + nf.) flowing into the nonlinear capacitor is Wan = Vin,almn® 
+ Va.n*Im.n, Where the asterisk denotes the complex conjugate. Then 
the general energy relations at the different frequencies are obtained. 
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These relations apply equally well to both the nonlinear capacitor and 
inductor, the two derivations being almost identical. 


2. Stability Theory and New Stability Criteria 


An excellent survey of recent work on theory of stability of solutions 
of nonlinear ordinary differential equations was given by W. Kaplan 
(101) in 1956. An extensive bibliography is appended, the publications 
being grouped under special topics, and over half of the references are to 
papers in Russian. Since there is no need of reminding the readers who 
know Russian about the Russian contributions, refs. 102, 103, and 104 
are given here as English translations are available. I. G. Malkin’s 
Methods of Lyapunov and Poincaré in the Theory of Nonlinear Oscilla- 
tions was published in Moscow, 1949; and his Stability of Motion was 
published in 1952.8 While Lyapunov (Liapounoff) first published his 
paper in 1892, a French translation by E. Davaux appeared in 1907, 
which was reprinted by the Princeton University Press in 1949. Re- 
prints of Lyapunov’s four papers in Russian appeared in 1950. A. I. 
Lur’s published his treatise Some Nonlinear Problems in the Theory of 
Automatic Regulation in 1951. An English translation is now available. 

According to Kaplan, for the work outside the Soviet Union, the 
following are important contributors: Bellman (105), Levinson (106), 
Wintner (107) in the United States, Miss Cartwright (108) in England, 
Haag (109) in France, Cesari (110) in Italy and the United States, 
Massera (111) in Uruguay. With regard to stability of forced vibra- 
tions, the recent contributions of refs. 112, 113, and 114 are worthy 
of note. 

Stability criteria in linear and nonlinear systems have been sug- 
gested from time to time. To cite the contributions at the Cranfield 
Conference 1951, L. S. Dzung (115) gave a stability criterion which is 
“more exact and more generalized” than the Nyquist criterion, R. W. 
Jones (116) applied the method of Liapounoff to give stability criteria 
for certain nonlinear systems, and A. Leonhard (117) proposed ‘relative 
damping”’ as criterion for stability and as an aid in finding the roots of 
a Hurwitz polynomial. (Some recent German contributions by K. 
Magnus, L. Cremer, H. Cremer, and F. H. Effertz are listed in refs. 
117b, c, d, e, and f.) 

For the stability criteria of nonlinear control systems, we shall 
cite two recent examples. K. Klotter and E. Pinney (118) gave a 
comprehensive stability criterion for forced vibrations in nonlinear 
systems as follows. Given the equation 


ax” + cx +c f(x’,x) = P cos ait. (50) 


8 There are four new Russian books called to the author's attention by Dr. B. Naumov. 
(See refs. 102b-—e). 
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Under very weak restrictions it is shown that the stability depends on 
the sign 0A/0@P, where A denotes the maximum value of x(t) within a 
period. The motion is stable if this derivative is positive ; it is unstable 
if it is negative. The plot of A, the maximum value of x(t) versus w,, 
the frequency of the driving term, or w,;? is known as a response curve. 
Once these response curves are established, one can proceed to study 
whether the motion corresponding to any one point on these curves is 
a stable one or not. Since Eq. 50 has no known exact solution, one 
often resorts to methods of integration which provide approximate 
results. Besides the perturbation and iteration methods, there are the 
averaging methods, of which the Ritz or Galerkin method is ay ex- 
ample (119). 

A new stability criterion for a class of nonlinear systems was more 
recently given by Ku and Wolf (120). According to the method of 
partition (94-96) the highest derivative x(t) of a k order nonlinear 
differential equation is expressible as a summation of w, ,t", as shown 
in Eq. 41. The response x(t) is then of the form 


= She. (51) 


n=0 


The Laplace transform of x(t) given by Eq. 51 is 
(52) 


The stability criterion is given as follows. The system is stable if and 
only if the convex singularity hull (CSH) of X(s) is entirely in the 
left-half plane. The CSH is constructed by joining the points of 
singularity of the transform X(s) in the s-plane. It may be noted that 
the singularities of X(s) are not necessarily poles. The degree of 
stability is also discussed in connection with the conjugate indicator 
function, which is called the Phragmen-Lindelof function. 


3. Optimum Nonlinear Control Systems 


According to Minorsky (121), some of the principal ways in which 
nonlinear relations may be used on control systems are: (1) to compen- 
sate for nonlinearities ; (2) to control the appearance or disappearance of 
an oscillation ; (3) for frequency de-multiplication ; (4) for synchroniza- 
tion; (5) for asynchronous excitation and quenching; and (6) in para- 
metric excitation. 

A theory of optimum nonlinear control (122), developed by R. 
Oldenburger, is devoted to systems with the differential equation 

_ P(D)e-Ta 


m' + f(m) = “O(D) g(c) — Lit) (53) 


7 
| 
4 
4 
{ 
is 


Feb., 1961.] THEORY OF NONLINEAR CONTROL 133 
where m denotes the deviation in the value of the controlled variable, 
c denotes the deviation in the value of the controlling variable, P(D) 
and Q(D) are polynomials in the derivative D = d/dt and with non- 
zero constant terms, e~74? denotes the dead-time operator, g(c) denotes 
a function of c, f(m) denotes a function of m, and L(t) is a given forcing 
function of time. Optimum transients for control systems can be 
obtained by nonlinear control functions depending on m. The results 
of theoretical studies have been verified experimentally with physical 
devices and automatically controlled systems. 

Design and analog-computer analysis of an optimum third-order 
nonlinear servomechanisms were reported by Doll and Stout (123). 
Nonlinear integral compensation of a velocity-lag servomechanism with 
backlash was made by Shen, Miller, and Nichols (124). Quasi-optimiza- 
tion of relay servos by use of discontinuous damping (125) was also 
reported. In the field of discontinuous automatic control, I. Fliigge- 
Lotz’s treatise (126) was published in 1953. She and Wunch (126a) 
showed that a system with discontinuous damping and spring coeffi- 
cients would transfer inputs well, if the change of the coefficients from 
one constant value to another was controlled by the deviation between 
input and output. Optimization by the phase-plane technique was 
given in refs. 70, 71, 72, and 73. Optimum response of discontinuous 
feedback control systems was analyzed by Nesline (127). A simple non- 
linearized control system was given by Zaborszky (128). Cosgriff and 
Emerling suggested optimizing control systems (129) using logic circuits. 

Compensation of feedback-control systems subject to saturation has 
been extensively studied by Newton (130). Nonlinear compensating 
networks were recommended by Mishkin and Truxal (131). A non- 
linear control system for wide range input signals was reported by Tou 
and Ku (132). Long and Ham (133) suggested conditional feedback as 
a new approach to linear and nonlinear feedback control. Compen- 
sating transient saturation by excess error storage was analyzed by 
Chang and Archbald (134). 


4. Sampled-Data Systems and Systems with Random Inputs 


The vast amount of literature on sampled-data systems cannot be 
adequately reviewed here (135-7). It can be said, however, that non- 
linearities in sampled-data control systems are now studied, and the 
theory of nonlinear control may well influence the design of sampled- 
data systems. As examples we may mention the recent papers by 
Kalman (138), Nease (139), Tou (140), and Mullin (141). 

The analysis of sampled-data systems owes much to the pioneer 
work of Ragazzini and Zadeh (142). Analysis of sampled-data servo- 
mechanisms has been done on a digital computer by Tostanoski (143) 
and on an analog computer by Chestnut (144) and others. It may be of 
interest to note that Boxer and Thaler (145) made use of the z-transform 
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in solving linear and nonlinear systems. Tsypkin’s Theory of Inter- 
mittent Regulation has been translated into English (146). 

As an introduction to the study of nonlinear systems with random 
inputs, we like to mention the papers by R. C. Booton, Jr. (54, 147), 
J. F. Barrett and Coales (148), W. E. Thomson (148b), J. C. West and 
Nikiforuk (149), K. B. Tuttle (150), D. W. C. Shen (151), K. Chuang 
and Kazda (151a). In ref. 54, Booton considered an amplitude-sensi- 
tive nonlinear element with a Gaussian input and derived a quasi- 
linearization representation of its response. ‘‘When the nonlinearity 
occurs as part of a feedback system, application of the quasi-lineariza- 
tion method is complicated by the fact that this method assumes 
knowledge of the probability density function of the input to the non- 
linearity. Since, in general, nonlinear operations on an input function 
change the shape of its probability density function, the inputs. to the 
nonlinearity and to the overall system do not have the same form of 
probability density function.”” Even if the input to the nonlinear 
element is assumed Gaussian, its rms value is unknown. Also the 
equivalent gain of the element is another unknown parameter, and two 
relations must be found to determine the equivalent gain. One of these 
relations is furnished by the quasi-linearization method. The second 
relation comes from analysis of the over-all quasi-linear system with 
the equivalent gain treated as a parameter. It is to be noted that the 
quasi-linear representation of the response is interpreted in terms of 
statistical characteristics of the actual response, which is described by 
the autocorrelation of the response and by the cross-correlation between 
the input and the response. 

Wiener’s theory (7) has been studied by Singleton (152), Ikehara 
(152b), Weinberg and Kraft (152c). The recent work of Bose (47, 
48) already mentioned in Sec. II-3 of this paper, is a further develop- 
ment of Wiener’s theory of nonlinear systems. Zadeh’s contribution 
to the theory of nonlinear systems (58, 153) is a parallel development 
and deals with optimum nonlinear filters. In ref. 58, it is shown that 
a general nonlinear two-pole of class N; is completely characterized 
by its responses to a family of step functions with amplitudes ranging 
over all real values. Several possible modes of realization of two-poles 
of class N,; are indicated. It is noted that, in order to determine the 
optimum filter within the Class N,, it is necessary to know the Second 
order probability density functions for the signal and noise. Optimiza- 
tion of nonlinear filters is given in ref. 153, and the extension of the theory 
to nonlinear multipoles is given in ref. 153a. 

Zadeh remarked rather optimistically (58) : 
The theory of nonlinear systems has made some progress recently toward the formula- 
tion of general theorems bearing on both the internal and external behavior of broad 


classes of nonlinear systems. 
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In making the above remark, he had in mind refs. 97, 98, 152a, 154, and 
154a. We might add ref. 155 to this group of papers. 

In connection with the optimum nonlinear filter, Zadeh (153b) 
pointed out two approaches: (1) one approach parallels Wiener’s ap- 
proach (to linear systems) by choosing the form or class of filters and 
then finding the optimum number of this class by minimizing the mean 
square error between the desired output and the actual system output ; 
(2) the other approach formulates an appropriate statistical criterion 
and then determines the optimum filter for this criterion with little or 
no restrictions upon the form of the filter. Both of these approaches 
yield equations for optimum filters in terms of higher-order statistics 
(distribution functions and correlation functions) of the input and the 
desired output. 

The new approach uses a polynomial functional method of charac- 
terization for the analysis and synthesis of nonlinear systems excited 
by white noise. This is the method of Wiener (7), Bose (48), and Bar- 
rett (148). This approach assumes that we have at our disposal an 
ensemble member of the input time function x(t) and the corresponding 
ensemble member of the desired output z(t). When these ensembles are 
available, the new approach makes measurements on them that yield 
optimum nonlinear systems (filters) directly, instead of first measuring 
the distributions and then solving design equations in terms of these 
measured values. 

Thus we see in the study of nonlinear systems with random inputs, 
we have started with information theory and correlation functions and 
now reached a method using ensembles and polynomial functions. 
We shall mention also that an Ensemble-Transform method has been 
recently developed by Ku and Wolf (156) for the study of nonlinear 
systems with random inputs (157). This method can be used with 
Taylor-Cauchy transforms for nonlinear systems or Laurent-Cauchy 
transforms for systems that can be described by differential-difference 
and sum equations. 


VII. CONCLUDING REMARKS 


From the above sections we see that the theory of nonlinear control 
systems has been developed rapidly in the recent years. It is the 
author’s hope that by reviewing the recent advances and interchanging 
new ideas at the first International Congress of IFAC on Automatic 
Control further progress can be promoted. 
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READER OPINION 
ON ELECTRIFICATION OF CARBON PLACK 


I have read the article on ‘Electrification of Carbon Black by 
Contact with a Metal Surface,’’ by Brasefield, which appeared in the 
October 1960 issue of the JOURNAL. Not enough work has been done 
in this field and it was a pleasure to learn of these studies. 

The electrical charge on cigarette smoke is a somewhat elusive thing 
to prove and I am curious to know more of this technique. I find it 
difficult to reach a conclusion regarding the electric charge especially 
if the only measurements made are those described in the article. 
There are questions on aging of the smoke particles and distinguishing 
between charges on smoke and gaseous ions. 

In work performed at Stanford University in the Electrical Engi- 
neering Department an effort was made to perfect the use of a radio- 
active probe in a Faraday cage. The results were quite disappointing. 
Several probes were investigated including polonium, tritium and 
radium D. Part of the difficulty can be attributed to the great vari- 
ability of the contaminates in the air and their distance from the radia- 
tion source in the space being tested. The polonium, alpha source, 
gave different results from the hard beta of radium D and different 
again from the soft beta of tritium. All three cause ionization of air, 
but because air is such a variable mixture, the space charge that is 
measured is greatly altered by the difference in ion mobility of charged 
submicroscopic particles. 

In my work with cigarette smoke I came to an opposite view about 
the electric charge after observing both the effect reported in the 
article and the direct effect of collecting smoke particles. The probe 
indicated a negative charge more often than positive but not con- 
sistently. Is it possible that the small positive ions become attached 
to cigarette smoke and the remaining small negative ions are measured 
by the collecting probe? 

Given time to age, the smoke particles acquire a positive charge 
from the small positive ions—leaving predominantly small negative 
ions with a high mobility. Since small ions have a mobility more than 
a thousand times greater than charged cigarette smoke, the charge 
sensed by the probe connected to the electrometer is negative. 

In separate experiments using cigarette smoke and an ion collector 
designed to measure ions of different mobility it was observed with fresh 
smoke taken directly through the cigarette, as would occur when a 
cigarette is smoked by a person, the smoke particles were charged 
equally positive and negative. The small ions in this set of experiments 
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were predominantly positive—assumed to be CO, and the few negatives, 
O,. Inasmall chamber where the smoke could age for a few seconds the 
smoke acquired a predominantly positive charge. The ion collector in 
these experiments was built with two stages and all measurements were 
made by recording positive and negative ion mobilities simultaneously. 

Differences were observed between smoke drawn through the ciga- 
rette and smoke taken from the lighted end. Also, smoke taken into 
the mouth and exhaled was relatively free of charge—showing a high 
adsorption of charge in the mouth and nasal passages. 

In summary, I question whether a Faraday cage test is valid for 
determining the electrical charge on cigarette smoke. A new instru- 
ment being developed at Stanford University may be the solution to the 
conflicting reports on particle charge and space charge. This instru- 
ment is designed to take the complete mobility spectrum of ions in air 
and automatically record the results. It offers an opportunity to 
identify some airborne submicroscopic contaminates by the resulting 
spectrum. A report was made on this new device before the Air Pollu- 
tion Symposium in New Orleans last March. However, the work is 
not yet complete and there are technical problems to be solved. 

J. C. BECKETT 
Electrical Engineer 
Palo Alto, Calif. 


Author's Reply. 


In reply to the above comments, may I first call attention to a report 
by B. Vonnegut and C. B. Moore entitled ‘“‘A Study of Techniques for 
Measuring the Concentration of Space Charge in the Lower Atmos- 
phere,’ published 31 January 1958 by Arthur D. Little, Inc., Cam- 
bridge, Mass., for the Geophysics Research Directorate of the Air Force 
Cambridge Research Center. The first part of this report, and par- 
ticularly page 48, presents in considerable detail the theory and use of a 
Faraday cage for measurement of space charge. It is possible that this 
presentation will clarify some of the questions raised. 

If the radioactive probe in the Faraday cage is at a positive potential 
with respect to its environment in a cloud of negatively charged smoke, 
negative ions produced by the alpha particles will be collected by the 
probe until its potential is lowered to that of its environment. The 
alpha particles, of course, produce as many positive ions as negatives, 
and these positive ions will effectively neutralize part of the negatively 
charged smoke. The effect of these positive ions on the initial elec- 
trometer reading will be small if the capacity of the electrometer is small. 
It is not surprising that results were different depending on whether 
the radioactive probe was an alpha emitter or a beta emitter, since 
the space potential measured by the probe is more localized in the 
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former case. However, the space charge should not be appreciably 
altered by the various probes if the electrometer capacity is small. 

We observed, as did Mr. Beckett, that when cigarette smoke is in- 
haled and held in the lungs for a number of seconds, the charge seems 
to dissipate because when the smoke is then exhaled into the Faraday 
cage, the electrometer shows only a small deflection. If a mouthful of 
smoke is discharged immediately into the Faraday cage, the electrom- 
eter generally shows a large negative deflection. However, if air is blown 
through the cigarette in such a way that the smoke from the lighted end 
drifts into the Faraday cage, we found that the smoke was occasionally 
positively charged. 

We performed another interesting experiment which I shall describe 
briefly. When it was first observed that cigarette smoke inside the 
Faraday cage was charged, the question arose as to whether the charge 
separation took place at the cigarette or at the walls of the Faraday 
cage as the smoke passed into the cage. We therefore constructed 
a second cage of the same dimensions as the first, but the walls of the 
second cage were made of 2-in. chicken wire instead of window screening. 
The cages were set up near each other and smoke was blown first into 
one and then into the other. Within the experimental error of these 
experiments, there was no difference in the electrometer readings of the 
two cages, suggesting that the smoke was already charged as it ap- 
proached the cage. 


The Faraday cage with polonium probe may be open to criticism 
as a quantitative instrument for measuring space charge on cigarette 
smoke but I can see no objection to its use as a qualitative instrument 
in experiments such as those we have been discussing. Nevertheless, 
I will be interested in learning more about the instrument under 
development at Stanford University. 


CHARLES J. BRASEFIELD 
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BOOK REVIEW 


OuTER SPACE PHOTOGRAPHY FOR THE AMA- 
TEUR, by Henry E. Paul. 124 pages, 
illustrations, 4% X8 in. New York, 
Amphoto, 1960. Price, $2.50. 


With the dawn of the space age there also 
was evident a rebirth of interest in celestial 
objects and astronomical phenomena. This 
interest was stimulated by the attempts of 
many lay people to photograph the Sputniks 
first and then the American satellites that 
swarmed into the sky. Asa result of this new 
development, astronomers in general and the 
staff of The Franklin Institute Astronomy 
Department, in particular, were deluged 
with calls requesting information on cameras, 
films, focal lengths, exposure times and other 
photographic terms in addition to the basic 
questions in astronomy. Thus, the reviewer 
and his colleagues found themselves in the 
unenviable position of trying to answer and 


interpret technical questions to a lay audience. 
It could be done but it took considerable time 


which was unavailable. It was at that time 
that we wished for an easily understandable 
book on the photography of celestial objects. 

In Outer Space Photography for the Amateur, 
my very good friend Dr. Henry E. Paul has 
written a book with authority and balance. 
In an inexpensive, small book, the author has 
explored the entire subject of celestial pho- 
tography. It is written in a lucid, non- 
technical fashion to permit easy reference and 
provide the essential information to permit 
the photographing of any object in the sky. 

As an example of the content, we find that 
for satellites we need a fast lens F :4.5 or faster 
with a 40 to 60 degree field and the very 


fastest film. You just open the shutter and 
record the trace of the satellite. 

Do you want to photograph the moon? 
You need a long focal length lens or the image 
is too small to use. Also, you need a yellow 
filter to correct for the color-error in the lens. 
Dr. Paul recommends as a starting point you 
photograph the quarter moon at F :8 at 1/50th 
of a second on Panatonic-X film. 

The sun is a more difficult subject and re- 
quires extensive experimentation with a series 
of filters. Using a 100,000 filter, the author 
recommends 1/60th of a second at F:11 on 
the slowest film available. 

As far as stars are concerned, you can get 
snapshots of these. Even color photographs 
are possible. What about a drive for the 
telescope? The author shows that the motion 
of the star image in one second of time is 
0.0007X focal length in inches. Comets, 
Meteors, Nebulae, Aurorae are all covered 
in this work. 

The author is a recognized master in the 
design, manufacture and use of optical 
equipment. His versatility and knowledge 
have been distilled and are contained between 
the covers of his superlative new book. I 
will personally recommend this book to the 
many questioners who desire the answers to 
problems in astronomical photography. At 
twice the modest cost of $2.50 this would bea 
striking bargain in books. At this price, no 
one with a telescope or camera can afford to 
be without it. 

I. M. Levitt 
The Fels Planetarium 
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BOOK NOTES 


THE CONTROL OF MULTIVARIABLE SYSTEMS, 
by Mihajlo D. Mesarovié. 112 pages, 
diagrams, 6 X9 in. New York, John 
Wiley & Sons, Inc.; Cambridge, The 
Technology Press; 1960. Price, $3.50. 


This significant contribution to the theory 
of control systems is divided into two parts: 
Part I covers interrelations of multivariable 
systems and Part II is on the synthesis of 
such systems. The author’s thesis is that 
control systems should have a multivariable 
foundation rather than a single-variable one. 
He developed this theory in 1958-1959 while 
he was a Sloan Foreign Post Doctoral Fellow 
at MIT. He uses a V-canonical structure 
with a feedback controller as the model for 
his multivariable theory and outlines situa- 
tions in which combinations of the P, V, 
and H structures are suitable. Readers of 
the JouRNAL will recall Dr. Mesarovié's 
paper in the April 1960 issue, on “Dynamic 
Response of Large Complex Numbers.” A 
point of interest to the JoURNAL is that the 
list of references includes 33 items, 4 of which 
are JOURNAL papers, 15 are books, 5 are 


reports or unpublished papers, 8 are from 
regularly published society transactions, and 
1 is from a Russian journal. 


SPACE TRAJECTORIES, edited by the Technical 
Staff of the Research Division of Radiation 
Incorporated. 298 pages, diagrams, 6 X 9 
in. New York, Academic Press Inc., 
1960. Price, $12.00. 


This volume is the published proceedings 
of a two-day symposium held in December 
1959, to survey current space trajectory 
analysis. Sponsors of the symposium were 
the Advanced Research Projects Agency, 
The American Astronautical Society, and 
the Research Division of Radiation Incor- 
porated. The fifteen papers comprising the 
volume deal with military space, satellite 
and probe orbit computation, trajectory 
constants, trajectory computation, computer 
programs for space missions, space-borne 
computer design, re-entry trajectories, and 
space-trajectory instrumentation. The fif- 
teenth paper on computer programming 
methods is a composite of six different pro- 
grams. Discussion follows some of the papers; 


some have up-to-date lists of references. 
This is a valuable reference work for the 
fast-growing science of space exploration. 


CONTRIBUTIONS TO THE THEORY OF NoN- 
LINEAR OSCILLATIONS, edited by L. Cesari, 
J. LaSalle and S. Lefschetz. 286 pages, 
diagrams, 7 X 10 in. Princeton, Princeton 
University Press, 1960. Price, $5.00 
(paper). 

Number 45 in the Princeton Annals of 
Mathematics Studies, this is Vol. V of the 
“Contributions to the Theory of Nonlinear 
Oscillations.” The variety of topics covered 
in Vol. V includes optimal control, Poincaré’s 
Perturbation Method, fixed point theorem, 
non-descriptive algebra, Lagrange stable 
motions, asymptotic stability in 3-space, 
and periodic solutions of various systems. 
Twelve of the thirteen papers are in English; 
the other is in French (on a system of tra- 
jectories of dynamic systems). 


INTRODUCTION TO QUANTUM MECHANICs, by 
R. H. Dicke and J. P. Wittke. 369 pages, 
diagrams, 6 X9 in. Reading (Mass.), 
Addison-Wesley Publishing Co., Inc., 1960. 
Price, $8.75. 


This text is designed for an introductory 
course in nonrelativistic quantum mechanics 
at the graduate level, but it may be used for 
advanced undergraduates. Prerequisites are 
a knowledge of calculus and differential 
equations and a course of basic physics 
including classical mechanics. The final 
chapter deals with quantum-statistical me- 
chanics, a subject previously neglected in 
texts on quantum theory. 


MovecuLaR DistILLaTion, by G. Burrows. 
211 pages, diagrams, 54 X 8} in. New 
York, Oxford University Press, 1960. 
Price, $5.60. 


This book is primarily intended for the 
laboratory research worker. The underlying 
principles are dealt with according to simple 
kinetic theory and also by means of the 
standard terminology used for distillation 
processes. Engineering problems associated 
with the production and maintenance of the 
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necessary degree of vacuum are given atten- 
Gas and vapour flow at low pressures 
Examples 


tion. 
also receive detailed treatment. 
of degassers and molecular stills are described 
together with typical applications. 


SuRFACE by S. Tolansky. 
296 pages, illustrations, 6 X 8} in. New 
York, Interscience Publishers, Inc., 1960. 
Price, $9.00. 


This work aims to illustrate the wide field 
open for study of this subject and to present 
a good deal of original matter. It wishes to 
show what can be achieved by precision 
multiple-beam interferometry when applied 
to the microtopographies of surfaces. It is 
a full presentation of the theory and practice 
of those interferometric techniques that the 
author himself has developed. Some 350 new 
interferograms are given, which demonstrate 
a powerful three dimensional microscopy; 
in the up-down direction magnifications of 
500,000 and resolutions of even quite small 
crystal lattice spacings have been achieved, 
yet by very simple means. The book’s many 
topics include microtopographical studies 
on both natural and cleavage faces of a 
variety of crystals; the etching of crystals; 
studies of oscillating quartz crystals; the 
finish of machined and polished surfaces; the 
study of hardness test indentations; erosion 
effects produced by high speed impact; and 
erosion due to electric sparks. It should be of 
special interest to industrial and academic 
physicists, crystallographers, mineralogists, 
engineers, chemists and metallurgists, to 
electro-technics specialists and to the dia- 
mond, machine tool and watch industries. 


POLYPROPYLENE, by Theodore O. J. Kresser. 
268 pages, illustrations, 5 X 7} in. New 
York, Reinhold Publishing Corp., 1960. 
Price, $6.50. 


All phases of polypropylene properties, 
and applications are brought 
It offers 
information on this increasingly important 
thermoplastic at a level where it can be readily 
understood by readers without special training 


production 
together in this concise handbook. 


in polymer chemistry. Particular emphasis 
has been placed throughout the book on the 
relationship between the development of 
polypropylene and the future of high polymers 


Book 


in general. The author has also included 
sufficient information on the basic principles 
of stereochemistry, to point up the unique 
qualities of polypropylene. Included as a 
special feature is a chapter on chemistry 
that describes the work of Professor Guillo 
Natta on the discovery and development of 
polypropylene. 


PRINCIPLES OF FEEDBACK CONTROL, by 
Charles H. Wilts. 271 pages, diagrams, 
6X9 in. Reading (Mass.), Addison- 
Wesley Publishing Co., Inc., 1960. Price, 
$8.75. 


This treatment of the analytical methods 
used in the design of feedback systems is 
intended for use as a textbook in courses 
at the advanced undergraduate-graduate 
level, and for practicing engineers who wish 
to study the principles of feedback control. 
The reader is presumed to have a good back- 
ground in ordinary differential equations 
and complex variable theory. Although 
desirable, the student is not required to have 
any formal training in the use of the Laplace 
transform for solving linear differential 
equations. The main theme of the book is the 
stability problem. Attention is focused on the 
fundamental aspects of the problem, to the 
Routh-Hurwitz and Nyquist tests, to the 
generalized concept of a root locus, and finally 
to the problem of creating a satisfactory 
degree of relative stability in systems that are 
either unstable or insufficiently stable. 


On THERMONUCLEAR WAR, by Herman Kahn. 
651 pages, 6 X 9} in. Princeton, Princeton 
University Press, 1960. Price, $10.00. 


The author, whose ‘‘briefings’’ on the hor- 
rors of total war have served to stir up mili- 
tary and civilian controversy, prepared this 
volume while he was a Visiting Research 
Associate at Princeton. In his preface, 
Mr. Kahn says “This book examines the 
military side of what may be the major 
problem that faces civilization, comparing 
some of the alternatives that seem available 
and some of the implications in these choices.”’ 
He has restricted the discussion to the 
deterrence and waging of all-out thermo- 
nuclear wars between the United States and 
the Soviet Union. He dedicates his book to 
“the goal of anticipating, avoiding and allevi- 
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ating crises.'’ There are twelve chapters 
and five appendices, covering such topics as 
genetic consequences, the SAC, problems of 
arms control, nonmilitary and _ military 
deterrents, Civil Defense programs, timing 
and coordination problems. In tables of data 
and through dispassionate analysis, Mr. 
Kahn makes it only too clear that unless 
we understand that the problem is one of sheer 
survival and are willing to plan for and spend 
enough money adequately to prepare for all 
eventualities, we cannot hope to avoid 
catastrophe. Mr. Kahn's main purpose in 
writing the book is to show that it is possible 
to decrease the probability of catastrophic 
nuclear war and to lessen its destructive 
effects if it does occur. The book should be 
read by all Congressmen, state officials, and 
those who are responsible for formulating our 
national defense policy. All readers will not 
agree with Mr. Kahn's analysis, but they 
should find it difficult to refute his objective 
arguments. 


INTRODUCTION TO MECHANICS, MATTER AND 
Waves, by Uno Ingard and William L. 
Kraushaar. 672 pages, diagrams, 6 X 9 in. 
Reading (Mass.), Addison-Wesley Pub- 
lishing Co., Inc., 1960. Price, $9.75. 


This elementary textbook, intended for 
use in the first year of a two-year course in 
general physics for students of science and 
engineering, is based on such a course taught 
at MIT by the authors. Basic concepts of 
mass, momentum, energy, angular momentum 
and their conservation laws are stressed. 
The authors’ aim has been to maintain a 
balance between experiment and inductive 
reasoning. Sample chapter titles are: Inertial 
Mass, Impulse and Work, Examples of Forces 
and Motion, Temperature and Heat, Atoms 
and Molecules, Some Properties of Matter, 
Mechanics of Fluids, and Transverse and 
Longitudinal Wave Pulses. Numerous prob- 
lems are included; answers are provided for 
the odd-numbered ones. 


Tue CHEMYSTRY AND TECHNOLOGY OF FER- 
TILIZERS, edited by Vincent Sauchelli. 
ACS Monograph 148. 692 pages, illustra- 
tions, 6X9 in. New York, Reinhold 


Publishing Corp., 1960. Price, $18.00. 


This monograph covers the major changes 
the chemistry and 


and developments in 
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technology of the fertilizer industry since 
1935. In 24 chapters contributed by experts 
from government, industry and agriculture, 
the book covers raw materials, manufacture, 
processing and use of chemical fertilizers. 
Special chapters are devoted to the caking 
of mixed fertilizers, materials-handling equip- 
ment and gaseous effluents. The editor, in 
his foreword, states that chemical fertilizers, 
although relatively a new agricultural devel- 
opment, “are one of the most strategic 
weapons of modern agriculture.” The 
monograph brings up-to-date the information 
on this strategic weapon and suggests possible 
progress for the future. 


Continuous GEOMETRY, by John von Neu- 
mann. 299 pages, 6 X9 in. Princeton, 
Princeton University Press, 1960. Price, 
$7.50. 


Continuous geometry was invented by John 
von Neumann in 1935. It is a remarkable 
extension of projective geometry to the non- 
finite dimensional case, just as Hilbert and 
Banach spaces are such an extension for 
vector spaces. Now published posthumously, 
these notes represent lectures given in 1935 
and 1937 at Princeton University, and at the 
Institute for Advanced Study. The text 
is divided into three parts. Part I develops 
the axioms of continuous geometry, dimen- 
sion theory and for irreducible case the dimen- 
sion function D (a). Part II presents von 
Neumann's discovery that every comple- 
mented modular lattice of order >4 can be 
represented as the lattice of all principal right 
ideals of a suitable ring. Part III begins the 
dimension theory for the reducible case of 
continuous geometry and gives the funda- 
mental decomposition techniques, the notion 
of central cover and the unrestricted transi- 
tivity of perspectivity. 


From DvuaLism TO UNITY IN QUANTUM 
Puysics, by Alfred Landé. 114 pages, 
54 X 84 in. New York, Cambridge Uni- 
versity Press, 1960. Price, $3.75. 


The author, in challenging the view that 
the principles of wave-particle duality, 
complementarity and uncertainty are funda- 
mental to micro-mechanics, sets out to reduce 
the dualistic quantum phenomena to still 
more elementary general principles of sym- 


; 
33 
2 
; 


152 


metry, invariance, continuity and the like. 
He thus develops the current Bohr-Heisen- 
berg account, resting on Born’s statistical 
interpretation, into a unitary theory. The 
book begins with a philosophical discussion 
of determinism and chance, in which the 
uncertainty of an individual event in a 
statistical group is linked to the postulate 
of macro-physical cause-effect continuity. 
Chapter Two deals with the states of a micro- 
physical object in general, and the transition 
probabilities from state to state in reaction 
to macro-physical measuring instruments. 
Chapter Three suggests that the quantum 
law of probability interference is not a mere 
oddity, but a consequence of the nature of 
probability itself as a lawful connection 
between events. Chapter Four gives a 
deduction of the most characteristic quantum 
feature—the wave-like relation between co- 
ordinates and momenta—from non-quantal 
theorems. The book concludes with a 
critical survey of present fashions in quantum 
philosophy and language. 


AN INTRODUCTION TO HOMOLOGICAL ALGEBRA, 
by D. G. Northcott. 282 pages, diagrams, 
6 X9 in. New York, Cambridge Univer- 
sity Press, 1960. Price, $8.00. 


This new mathematical subject has emerged 
during the past ten years. Because of its 
fundamental nature it is relevant to many 
branches of pure mathematics, including 
number theory, geometry, group theory and 
ring theory, and it is a subject which mathe- 
maticians wishing to keep informed about 
major developments, cannot neglect. The 
author aims to introduce homological ideas 
and methods, and to show some of the results 
which can be achieved. The early chapters 
provide the results needed to establish the 
theory of derived functors, and to introduce 
torsion and extension functors. The new 
concepts are then applied to the theory of 
global dimension, in an elucidation of the 
structure of commutative Noetherian rings 
of finite global dimension, and in an account 
of the homology and cohomology theories 
of monoids and groups. A final section is 
devoted to comments on the various chapters, 
supplementary notes, and suggestions for 
further reading. This work will be useful as a 
reference for specialists as well as a helpful 
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and lucid account for those about to begin 
research. 


DiGITaAL APPLICATIONS OF MAGNETIC DE- 
vicEs, edited by Albert J. Meyerhoff. 
604 pages, diagrams, 6 X 9 in. New York, 
John Wiley & Sons, Inc., 1960. Price, 
$14.00. 


The first comprehensive text on digital 
magnetic devices, this volume will be a boon 
to students and professional computer engi- 
neers. It is based on the research programs 
sponsored by Burroughs Corporation on 
magnetic-core devices and magnetic-shift 
registers. The book covers theory, electrical 
design and logical design necessary for the 
development of digital systems using the 
magnetic core as the basic circuit element. 
The text is divided into the following parts: 
Fundamentals, Parallel Magnetic Pulse Am- 
plifiers, Delay Parallel Magnetic Pulse 
Amplifiers, Series Magnetic Pulse Amplifiers, 
Memories, Transistor-Magnetic Pulse Am- 
plifiers, and Other Techniques. A good 
feature is the inclusion of a list of symbols at 
the beginning of each part. 


ANNUAL REPORTS ON THE PROGRESS OF 
CHEMISTRY FOR 1959. Vol. LVI, 476 pages, 
diagrams, 54 X 84in. London, The Chem- 
ical Society, 1960. Price, 2 pounds, post 
free. 


The 1959 volume of these Annual Reports 
is divided into five sections: General and 
Physical Chemistry (9 papers); Inorganic 
Chemistry (3 papers); Organic Chemistry 
(11 papers); Biological Chemistry (5 papers); 
and Analytical Chemistry (7 papers). The 
several papers fall into the same general 
classifications as last year’s volume, but they 
either cover different phases of the subject, 
or extend the coverage. One entirely new 
topic (for the Reports) occurs in the Bio- 
logical section—the chemistry of bone and 
tooth minerals. As in past volumes, the 
reference lists are quite valuable, covering 
for the most part only 1959 entries. Refer- 
ences are given as footnotes throughout the 
text ; an alphabetical index of authors referred 
to is included at the end of the volume. 


THE Story OF STEREO: 1881, by John Sunier. 


160 pages, illustrations, 5} X 8} in. New 
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York, Gernsback Library, Inc., 1960. 
Price, $5.00 (hard cover); $2.95 (soft 
cover). 


The history and development of stereo 
from its earliest beginnings in the last century 
right up to the present are given in this book. 
It is also designed to aid in an understanding 
of the principles behind stereophony and what 
the medium has accomplished. The applica- 
tions of stereo on film, tapes and discs, and 
also in broadcasting have been covered by 
the author. He concludes by describing 
stereo techniques in the home, in business 
and industry and in medicine. 


ELECTRONIC MAINTAINABILITY, edited by 
F. L. Ankenbrandt. Vol. 3. 312 pages, 
illustrations, 6 X 9 in. Elizabeth (N. J.), 
Engineering Publishers, 1960. Price, $10.00. 


This book is based on the Third Electronic 
Industries Association Conference on Main- 
tainability of Electronic Equipment, held in 
San Antonio, Texas. Twenty-nine papers 


are herein presented giving an authoritative 
discussion of the developments which have 
been changing maintainability from the 


vaguely practiced art of yesterday to the 
rigorous science of today and tomorrow. 
Subject matter ranges from maintenance 
problems of space flight to the determination 
of adequate working space for electronic 
technicians. The contributors are among 
the foremost experts in the field. Although 
concerned primarily with military electronic 
equipment, much of the material is equally 
applicable to more prosaic products such as 
automobiles, household appliances, plumbing 
fixtures, and what have you. 


Orr-THE-RoaD Locomotion, by M. G. 
Bekker. 220 pages, illustrations, 6 X 9 in. 
Ann Arbor, University of Michigan Press, 
1960. Price, $10.00. 


Mr. Bekker, whose interesting paper in 
this JouURNAL for May 1957 introduced many 
of our readers to his work in off-the-road 
locomotion, presents a comprehensive treat- 
ment of this new approach to an age-old 
transportation problem. Mr. Bekker’s main 
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contribution to terramechanics is to relate the 
vehicle as a whole to the terrain it must 
operate on. Methods of predicting vehicle 
performance, soil and snow mechanics, tracks, 
wheels and tires are all treated in this com- 
pact volume, with emphasis on the fact that 
sound theory must be used in designing 
equipment in these days when “‘cut-and-try” 
methods are too slow and expensive. A 
wealth of graphs and photographs is most 
helpful, especially when the author delves 
into material not familiar to the reader. 
In the chapter describing the effective spaced 
link track, however, the eight photographs 
of Fig. 75 have been reduced so greatly in 
reproduction that all detail has been lost in a 
meaningless blur. 


CLASSICAL ELECTRICITY AND MAGNETISM, by 
E. S. Shire. 396 pages, diagrams, 54 X 8} 
in. New York, Cambridge University 
Press, 1960. Price, $7.50. 


This volume based on classical electricity 
and magnetism bridges the gap between the 
elementary account as taught in schools and 
a mathematical treatment based on Maxwell’s 
equations. It deals with the “classical,” 
that is, macroscopic, aspects of the subject 
to the standard of a general degree course 
in physics or engineering or of the first part 
of an honors degree course. The author 
combines an experimental with an historical 
approach. There are chapters on the meas- 
urement of electrical quantities and on applied 
electricity, the latter including introductory 
accounts of electronics and of machines for 
accelerating charged particles. Simple vector 
methods are used but vector calculus equiva- 
lents of the more important equations are 
confined to footnotes and appendices. An- 
swers are given to the wide range of examples. 
In the belief that a physicist should be able 
to adopt whichever system of units is most 
convenient for the purpose in hand, the author 
has not tied his book to any one system, but 
has made it equally applicable to ‘‘rational- 
ized” and m.k.s. and c.g.s. 
systems. The symbols used are those recom- 
mended by the International Union of Pure 
and Applied Physics. 
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PUBLICATIONS RECEIVED 


LAPLACE TRANSFORMATION, by William T. Thomson. Second edition, 255 pages, diagrams, 
6 X 9in. Englewood Cliffs, Prentice-Hall, Inc., 1960. Price: $10.00 (trade) ; $7.50 (text). 


TuBEs AND Circuits, by George J. Christ. 192 pages, diagrams, 54 X 84 in. New York, 
Gernsback Library, Inc., 1960. Price: $5.00 (hard cover); $3.45 (paper). 


MECHANICAL-ELECTRICAL EQUIPMENT HANDBOOK FOR SCHOOL BuILDINGs, by Harry Terry. 
412 pages, illustrations, 6 X 9 in. New York, John Wiley & Sons, Inc., 1960. Price, 


$9.50. 


ACTIVATION ANALYSIS HANDBOOK, edited by R.C. Koch. 219 pages, 8 X 10}in. New York, 
Academic Press Inc., 1960. Price, $8.00. 


FuNDs FOR RESEARCH AND DEVELOPMENT IN INDUSTRY, 1957. 119 pages, diagrams, 8 X 10 
in. National Science Foundation Survey 60-49. Washington, Government Printing 
Office, 1960. Price, $0.65 (paper). 

MATHEMATICAL TABLES. VOLUME 3. TABLES OF GENERALIZED EXPONENTIAL INTEGRALS, 
by G. F. Miller. 43 pages, 8$ X 10}in. London, Her Majesty's Stationery Office, 1960. 
Price, 7s, 6d. American agents, British Information Services in New York; price, $1.43 
postpaid. 

SCIENTIFIC AND TECHNICAL PERSONNEL IN AMERICAN INDUSTRY. Report on a 1959 Survey. 
66 pages, 8 X 10}in. National Science Foundation Survey 60-62. Washington, Govern- 
ment Printing Office, 1960. Price, $0.45 (paper). 

Les STRUCTURES DE COMMUTATION A m VALEFURS ET LES CALCULATRICES NUMERIQUES, by 
Midhat J. Gazalé. 78 pages, 6} X 94 in. Paris, Gauthier-Villars, 1959. Price, $3.04 
(paper). 
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NATIONAL BUREAU OF STANDARDS NEWS 
THE AMOS IV COMPUTER 


The National Bureau of Standards in cooperation with the U. S. 
Weather Bureau has developed a specialized digital computer! for the 
Weather Bureau to use as a research tool in exploring the concept of 
the automatic weather station. The AMOS IV computer receives 
data from weather-sensing instruments and processes these data through 
such functions as sampling, comparing, selecting a maximum, and arith- 
metic operations. The results are transmitted via teletype to a central 
forecasting station and to other airport weather stations. Values of two 
quantities recently developed as aids to air safety—runway visual 
range and approach light contact height—are given by the machine 
through automatic table look-up. 

For a number of years, the Weather Bureau has been appraising the 
possibilities of an automatic weather station. Such stations could be 
widely distributed, and would be especially useful in relatively inacces- 
sible locations that are important sources of early data on meteorological 
activity. The various developmental prototypes of this concept have 
been called AMOS (Automatic Meteorological Observation Station) ; 
the current version, containing transistorized packages, is AMOS IV. 
This model was designed and built by Paul Meissner and J. A. Cunning- 
ham of the NBS data processing systems laboratory and by C. A. 
Kettering of the U. S. Weather Bureau. It is an outgrowth of previous 
work done by NBS for the Weather Bureau that resulted in a special 
computer? for processing cloud-height signals from a ceilometer. The 
ceilometer was intended for use with the AMOS III. 

Several of the input quantities to the AMOS computers, such as 
cloud height and precipitation, cannot be satisfactorily represented by 
instantaneous values but must be time-averaged. Varying amounts 
of data processing must therefore be associated with the different 
instruments measuring these quantities. In the AMOS III concept, 
several complex units were required for these functions. Although 
many of the functions were similar, the hardware was not minimized 
because of a diversity of design that resulted from the isolated develop- 
ment of the individual units. Analysis of the over-all system indicated 
that a considerable reduction could be made in hardware and therefore 
in maintenance. 

In AMOS IV, the automatic weather station is built around a single 
small, general-purpose computer designed especially for this application. 


1“A Computer for Weather Data Acquisition,” by P. Meissner, J. Cunningham, and C. 
Kettering (to be published Proc. of EJCP, 1960). 
2 “Cloud-Height Data Analyzer,” NBS Tech. News Bull., Vol. 43, p. 180 (1959). 
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The computer receives data from the input instruments at any desired 
interval. These data are suitably processed and arranged in a specified 
order for teletype transmission in a variety of message formats and at 
various speeds. The computer also operates local and remote displays. 
Much latitude is available for research into the most desirable form of 
data processing because of the inherent flexibility of the internally 
programmed machine. 

The computer continuously monitors new input data while simul- 
taneously processing data already entered and transmitting messages 
on command. Among the input quantities which the AMOS IV 
computer can handle are temperature, dew point, wind speed and direc- 
tion, atmospheric pressure, precipitation, transmissivity, and cloud 
height. Input data can be received directly from the instruments in 
the simplest possible form, such as analog voltage, current, or resistance ; 
and pulse rate or contact closure. Information may also be received in 
coded form, such as the Gray binary code frequently used with shaft- 
position encoders. The nature of the weather instruments and of the 
quantities measured limits the input data to 2 or 3 decimal digits for 
the most part; word size is therefore 3 digits plus sign. Double-pre- 
cision methods are available for those few instances requiring greater 
accuracy. Communication with the machine is via electric typewriter 
or punched tape. 

The computer circuitry is based on transistorized plug-in assemblies* 
designed at NBS for a variety of data-processing applications. These 
50-kc packages perform flip-flop, analog switch, and gating circuitry 
functions, as well as others. 

To store data, the machine uses a magnetic drum operating at 1800 
rpm that carries 100 general storage channels of 100 words each and 
has space for 100 additional channels. Several dual-head channels are 
available for simultaneous read-in and read-out of incoming data, out- 
going messages, etc. The magnetic drum provides the extensive storage 
capacity required for the table look-up involved in the calculations of 
runway visual range and approach light contact height. About 35 
tables are stored on the drum ; each table has about 90 three-digit values. 

One set of these tables contains the data on runway visual range 
(RVR), i.e., the distance along the runway visible to a pilot from the 
point of touchdown—generally 1000 to 6500 feet, depending upon run- 
way illumination (natural and artificial) and atmospheric conditions. 
The primary input for the RVR determination is a transmissometer 
reading. The computer continuously monitors this reading and ‘‘looks 
up” the proper corresponding value of RVR, which is then displayed 
locally and inserted into the teletype message. 

The other set of tables contains the data on approach light contact 
height (ALCH), i.e., the height from which the pilot can identify the 


3 “Packaged Switching Circuits,”” NBS Tech. News Bull., Vol. 43, p. 184 (1959). 
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approach lights. ALCH is affected by background illumination level, 
atmospheric conditions, and the intensity of the approach lights, which 
are set in accordance with prevailing conditions. If limiting conditions 
are indicated by either low clouds, as shown by the ceilometer, or by 
fog or snow, as sensed by the transmissometer, a value of ALCH based 
on the interfering factor is obtained. If both factors are present, 
two calculations are made; the machine then determines and displays 
the lower value. Since there is a statistical uncertainty in this type of 
information, two values of altitude are presented. The higher altitude 
is that at which the pilot has a 20 percent probability of seeing the 
approach lights; the lower altitude is that at which the probability is 
90 percent. 


YOU CAN ADVANCE SCIENCE EDUCATION 


Today, more than ever before in its 135-year history, there is vital need for 
The Franklin Institute to effectively promote education in science and technology. 
It is imperative that we meet this challenge by providing adequate educational op- 
portunities in these fields. This requires vision, objective planning, and money. 


We have more than enough of the first two requisites, but far too little of the third. 

Our programs are aimed at professional scientists and industry, as well as the lay 
public and young people seeking inspiration and guidance in choosing a career. 
The Institute’s educational programs are impressive, for they begin with students in 
the early grades of our elementary schools and continue throughout an individual’s 
professional or industrial life. With more funds at our disposal, the scope and 
vigor of these activities could be greatly increased. 

The Franklin Institute is not richly endowed. It is a non-profit organization, 
depending for encouragement and support on an understanding public. Capable 
and conservative management assures wise administration of all funds. 

Your gift or bequest, large or small, will be deeply appreciated and will be used 
effectively to broaden the Institute’s educational usefulness. There is a warm 
satisfaction in giving financial support to an organization that has pioneered in, and 
is dedicated to, the advancement of science and technology through education. 

When property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

The Secretary of The Franklin Institute will gladly furnish you with additional 
information. Write to him at The Franklin Institute, Benjamin Franklin Parkway 
at Twentieth Street, Philadelphia 3, Pennsylvania. 


; 


THE FRANKLIN INSTITUTE 


MINUTES OF THE ANNUAL MEETING 
January 18, 1961 


The Annual Meeting of The Franklin Institute was held at 8:15 p.m. in the Lecture Hall. 
Wynn Laurence LePage, President, called the meeting to order with approximately 250 
members and guests in attendance. 

The President stated the minutes of the Stated Meeting of November 18, 1960 were 
published in the December issue of the JouURNAL. There being no additions or corrections, the 
minutes were approved as published. The minutes of the Stated Meeting of December 21, 
1960, he stated, will be published in the January issue of the JOURNAL and will be presented 
for approval at the next Stated Meeting. 

The President then called on the Chairman of the Tellers of Election, Dr. Joseph S. 
Hepburn, for the Report of the Tellers who had been appointed in accordance with the By Laws 
to count the votes for ten Managers to be elected. Dr. Hepburn reported that 2891 votes were 
cast, of which 2502 were legal and 389 defective, and the nominees received votes as follows: 


Morton Gibbons-Neff 2037 Wynn Laurence LePage 2325 
Rupen Eksergian 2192 Gaylord P. Harnwell 2332 
Francis J. Chesterman 2290 Emery W. Loomis 2254 
James Creese 2311 Henry M. Chance, I1 2263 


Henry B. Bryans 2253 Joseph A. Fisher 2264 
Kenneth Robert Atkins 1507 


The President then declared the following ten members were thereby elected Managers 
of The Franklin Institute for a term of three years: 


Morton Gibbons-Neff Wynn Laurence LePage 
Rupen Eksergian Gaylord P. Harnwell 
Francis J. Chesterman Emery W. Loomis 
James Creese Henry M. Chance, II 
Henry B. Bryans Joseph A. Fisher 


The President thanked Dr. Hepburn, Mr. Duerig and Mr. Neher for their services as 


Tellers of Election. 

Mr. LePage then called on the Secretary of The Institute to read the report submitted 
by the Trustees of the Elliott Cresson Medal Fund showing the financial transactions for the 
year, which, in accordance with the Deed of Gift that forms the basis of the Award, must be 
presented at the Annual Meeting of The Franklin Institute. The report follows: 


ELLIOTT CRESSON MEDAL FUND 


Report of the Trustees to The Franklin Institute 
for year ended December 31, 


PRINCIPAL 


Principal of the Fund (Representing the original $1,000 plus transfers thereto from eenanneaten Income 
and profit from disposition of investments) $4.3 


Invested as follows: 
$4,000 U. $4,023.79 


S. Treasury 24% bonds, due 6/15/1962-67 at book value 


Cash in Western Savings Fund Society 280.94 


INCOME 
Income of the Fund for 1960 
Interest on bonds (less amortization) 
Interest on bank balance 
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04.73 
$4,304.73 
$ 121.44 
9.14 $ 130.58 


Feb., 1961.1] New Books IN FRANKLIN INSTITUTE LIBRARY 159 


Expenses 
Medals, certificates, insurance and cases... . 232.89 


Expenses in excess of Income for the year 1960....... ‘ (102.31) 
Proceeds from Sale of $1,000 U. S. Treasury Bond. 936.25 


$ 833.94 


Unexpended Income at January 1, 1960. . 1,971.49 
Less: Value of bond sold..... 1,005.98 965.51 


Unexpended Income at December 31, 1960................. $1,799.45 


Inv as 
Treasury 24% bonds, due 6/15/1962-67 at book value... . 1,005.98 
Cash in Philadelphia Trust Company...... 793.47 


$i, 799.45 


s/ Henry B. Allen 
Henry B. ALLEN 


s/ Delmer S. Fahrney 
S. FAHRNEY 
Trustees 


January 18, 1961 


The President presented The President’s Letter, which will accompany the Annual Report 
of The Franklin Institute for 1960, which report will be distributed to the membership of The 
Institute early this year. The President's Letter contained highlights covering The Institute’s 
activities for 1960. 

Mr. LePage then presented Sir George Thomson, Master of Corpus Christi, Cambridge, 
and President of the British Association for 1960, as the speaker of the evening. Sir George 
spoke on “‘The Importance of Useless Science.’’ His talk on science as one of the great achieve- 
ments of the human intellect emphasized the debt our technology owes to free investigation 
and was illustrated with interesting slides. 

After thanking Sir George for his challenging talk, and before calling for questions, Mr. 
LePage presented Sir George with a copy, from our own Library, of his father’s book, titled, 
“The Electron in Chemistry’ which covers a series of five lectures delivered by Sir Joseph J. 
Thomson before The Institute in 1923. In this book, Mr. LePage said, he had taken the 
liberty of inscribing, ‘‘To Sir George Thomson in memory of your late father, Sir Joseph J. 
Thomson, whose scientific career was revolutionary in its successes,’’ and had signed as Presi- 
dent of The Franklin Institute. 

Following an interesting question period, the President adjourned the meeting at 9:45 P.M. 

WILLIAM F, JACKSON, JR. 
Secretary 


NEW BOOKS IN THE FRANKLIN INSTITUTE LIBRARY 


MATHEMATICS 


STATISTICAL THEORY OF COMMUNICATION, by Y. W. Lee. Wiley, 1960. 

AN INTRODUCTION TO THE THEORY OF NUMBERS, by I. M. Niven. Wiley, 1960. 

AN INTRODUCTION TO THE THEORY OF EXPERIMENTAL DesiIGNn, by D. J. Finney. University 
of Chicago Press, 1960. 

AUTOMATIC TRANSLATION, by D. I. Panov. 


Pergamon, 1960. 


ASTRONOMY 


Batchworth Press, 1959. 


LAROUSSE ENCYCLOPEDIA OF ASTRONOMY, by L. Rudaux. 

Moon Maps, by H. P. Wilkins. Macmillan, 1960. 

THE PLANETARY EQUATORIUM OF JAMSHID GuiyATH AL-Din AL-KAsui, with translation and 
commentary by E. S. Kennedy. Princeton University Press, 1960. 

INTRODUCTION TO THE MECHANICS OF THE SOLAR SysTEM, by R. Kurth. Pergamon, 1959. 

Source Book 1n Astronomy, edited by H. Shapley. Harvard University Press, 1960. 
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PHYSICS 


LA pes Gaz NeutREs Et Ionisés. Wiley, 1960. 
ENGINEERING THERMODYNAMICS, by N. A. Hall. Prentice-Hall, 1960. 
THEORETICAL ELASTICITY AND PLASTICITY FOR ENGINEERS, by D. E. R. Godfrey. Thames 
and Hudson, 1959. 
SEMICONDUCTORS, by R. A. Smith. Cambridge University Press, 1959. 
PLAsMA Puysics, by S. Chandrasekhar. University of Chicago Press, 1960. 
SPECIAL RELATIVviTy, by W. Rindler. Interscience, 1960. 
PROPERTIES OF MaTTER, by F.C. Champion. 3rd ed. Philosophical Library, 1959. 
LECTURES ON THEORETICAL RHEOLOGY, by M. Reiner. 3rd ed. Interscience, 1960. 
Low-TEMPERATURE TECHNIQUES, by F. Din. Interscience, 1960. 
CONFERENCE ON SILICON CARBIDE, Boston, 1959. Proceedings, edited by J. R. O’Connor. 
Pergamon, 1960. 
CHEMISTRY AND CHEMICAL ENGINEERING 


THE CHEMISTRY OF PROPELLANTS, edited by S. S. Penner. Pergamon, 1960. 

Boron FLuoRIDE AND Its CompounDs As CATALYSTS IN ORGANIC CHeEMisTRY, by A. V. 
Topchiev. Pergamon, 1959. 

IMPERFECTIONS IN CRysTALs, by H. G. van Bueren. Interscience, 1960. 

CHEMISTRY AND TECHNOLOGY OF FERTILIZERS, by V. Sauchelli. Reinhold, 1960. 

STABILIZATION OF FREE RADICALS AT Low TEMPERATURES, edited by A. M. Bass. U.S. Dept. 
of Commerce, National Bureau of Standards, 1960. 

SYMPOSIUM ON Spectroscopy. ASTM, 1960. 

Foop CHemistry, by L. Meyer. Reinhold, 1960. 

X-Ray PowpER PHOTOGRAPHY IN INORGANIC CHEMISTRY, by R. W. M. D’Eye. Academic 
Press, 1960. 

THEORY OF DETONATION, by I. B. Zel'dovich. Academic Press, 1960. 

GRowTH OF CRYSTALS, v. 2, edited by A. V. Shubnikov. Consultants Bureau, 1959. 

SELECTED READINGS IN GENERAL CHEMISTRY, by W. F. L. Kieffer. Division of Chemical 
Education of the American Chemical Society, 1959. 

MATERIALS SCIENCE 


MATERIALS IN NUCLEAR APPLICATIONS. ASTM, 1960. 

SYMPOSIUM ON RADIOISOTOPES IN METALS ANALYSIS AND TESTING. ASTM, 1960. 

SYMPOSIUM ON VISUAL AIDS FOR STANDARDIZING AND COMMUNICATING PRopuUCT APPEARANCE. 
ASTM, 1959. 

CORROSION PROBLEMS OF THE PETROLEUM INDUSTRY. Macmillan, 1960. 

ZONE REFINING AND ALLIED TECHNIQUES, by N. L. Parr. G. Newnes, 1960. 

SYMPOSIUM ON HypRAULic FLuips. ASTM, 1960. 


COMMITTEE ON SCIENCE AND THE ARTS 


(Abstract of Proceedings of Stated Meeting held Wednesday, January 11, 1961.) 
Hall of the Committee, 
Philadelphia, January 11, 1961. 
Mr, J. PutLip Evans in the Chair. 


The following report was presented for final action: 
No. 3325: The Bubble Chamber. 


This report recommended the award of an Elliott Cresson Medal to Donald A. Glaser, 
of Berkeley, California, ‘‘In consideration of his invention of the bubble chamber, a device in 
which the tracks of high energy ionizing particles and of the fragments of nuclear collisions can 
be photographed, thus providing a powerful tool for obtaining new information about ele- 


mentary particles and processes that give rise to them.” 
D. S. FAHRNEY, 


Secretary to Committee 
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Hypersonic Shock Tunnel.—Some 
of the conditions that will be en- 
countered by the first interplanetary 
space ship as it comes in for a landing 
on Mars or Venus are being simu- 
lated at the General Electric Research 
Laboratory in Schenectady, N. Y. 
Velocities many times the speed of 
sound, in atmospheres believed to be 
the same as those of Mars and Venus, 
are being produced in a shock tunnel. 
By studying the effects of such ve- 
locities on various shapes, scientists are 
getting a preliminary idea of what 
kind of space ship an astronaut 
will use in his first planet-hopping 
expedition. 

The General Electric shock tunnel 
has produced some of the highest 
flow velocities of air ever reported, 
up to 35 times the speed of sound, 
and is the first to provide two separate 
test chambers of identical high-speed 
air flow. The provision of two test 
chambers in a_ hypersonic shock 
tunnel permits many different types 
of measurements to be made under 
absolutely identical conditions for the 
first time. 

In conventional shock tunnels, only 
a limited number of measurements 
can be made during one test run, and 
it is extremely difficult to duplicate 
test runs exactly. 

The shock tube itself is simply a 
device for putting a large amount of 
energy into a quantity of gas. A 
mixture of hydrogen and oxygen, 
diluted with varying amounts of 
helium, is ignited, producing a con- 
trolled explosion. This “explosion” 


ruptures a stainless steel diaphragm 
separating the burning mixture from 
low pressure air in the rest of the tube. 
The shock thus produced travels the 
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100-ft. length of the tube, and 
bounces back from the far end. For 
a few thousandths of a second the 
air at this end of the tube is com- 
pressed to an extremely high pressure 
and temperature. Two nozzles open- 
ing off the end of the tube allow the 
highly compressed air to expand and 
flow over small models being tested. 
Depending on the size and shape of 
the nozzle opening and the strength 
of the initial shock, the flow of air 
over the test models may be as high 
as 35 times the speed of sound, or 
Mach 35. 

Many tests have been made at 
flows between Mach 20 and 25, with 
work now progressing in the range 
of Mach 25 to 35. Research at these 
speeds is of importance primarily 
for space vehicles, since Mach 35 is 
about the speed needed for escape 
from the earth’s gravitational field. 
Even high speeds, up to Mach 50 
in air, are probably attainable in the 
General Electric shock tunnel. In 
addition to high speed, practically 
any desired atmospheric conditions 
can be provided, and re-entry into any 
known atmosphere can be simulated. 

Shock tunnel studies provide some 
of the data necessary to predict how 
hot a space vehicle will get when re- 
entering the atmosphere, how the heat 
is distributed over the surface, and 
what pressures may be expected. 
Other types of shock tubes, such as 
the plasma jet or arc tunnel, are also 
being used by General Electric to 
study re-entry problems at its Mis- 
siles and Space Vehicle Department 
in Philadelphia. These studies are 
providing the basis for choice of 
materials and shapes for space vehicles 
of the future. 
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New Control System for Weather 
Satellite.—Tiros II, the nation’s new 
television-infrared weather ‘‘eye’’ sat- 
ellite, has achieved a major ‘“‘first”’ 
in space with a new and simple system 
that enables ground observers to shift 
its tilt in space by remote control 
for improved TV picture coverage of 
clouds around the Earth. 

The new orientation technique was 
developed by the Radio Corporation 
of America research and engineering 
team which designed and built both 
the Tiros I and Tiros II satellite-and- 
ground systems for the National 
Aeronautics and Space Administra- 
tion. It uses the effect of the Earth’s 
magnetic field in space to alter the 
attitude of the satellite upon com- 
mand, without the need for special 
propulsion devices such as rocket 
motors or compressed gas jets. Ap- 
plied for the first time in Tiros II, the 
system may have application to 
future satellite projects for stabiliza- 
tion and attitude control. 

The new system results from studies 
by RCA and government scientists 
of an unexpected gradual shift in the 
attitude of the first Tiros satellite 
under the influence of the magnetic 
field surrounding the Earth. In the 
first Tiros, which returned nearly 
23,000 useful TV cloud pictures to 
Earth following its launching last 
April 1, these magnetic forces caused 
the satellite to tilt gradually away 
from the predicted position of its 
axis in space. 

In Tiros II, the forces are being 
harnessed by a technique in which 
a controllable magnetic field is gen- 
erated around the satellite itself by a 
coil of wires around the lower sides of 
the vehicle. Interacting with the 
Earth’s magnetic field, this gives 


ground observers an invisible “hand” 
with which to tilt the satellite on 
command in order to obtain a more 
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advantageous angle for the sensors 
and the solar power supply. In 
developing the novel technique, RCA 
engineers devised a large spherical 
test “‘cage’’ in which the complete 
satellite was placed and rotated in 
a magnetic and solar environment 
similar to those in space. 

Equipped with its new orientation- 
control system and with newly-devel- 
oped infrared instruments to measure 
the emission and reflection of solar 
heat by the Earth and its atmosphere, 
Tiros II is the second step in the 
experimental weather satellite pro- 
gram being conducted by NASA to 
study the feasibility of regular satellite 
weather operations. 

The first step was Tiros I, whose 
three months of operation during last 
yéar shed new light on global weather 
patterns and demonstrated the utility 
of satellites for weather observation. 
Both the Tiros I and Tiros II sys- 
tems and their associated ground net- 
works of command and receiving sta- 
tions have been developed for NASA 
by RCA’s Astro-Electronics Division, 
Princeton, N. J. 

In the satellite itself, the major 
innovations of Tiros II are the 
orientation-control system and the 
infrared equipment to observe the 
heat balance of the Earth and its 
atmosphere. The infrared package, 
developed by NASA, includes a five- 
channel detector to measure selected 
portions of the infrared spectrum 
around the Earth. These provide 
information about the amount of 
reflected sunshine, the separate and 
total amounts of radiation from the 
Earth and its atmosphere, and, for 
reference, the amount of visible light 
reflected back into space. 

It has been pointed out that in 
spite of the addition of new features 
to Tiros II, its total weight has been 
kept down to about 280 lb. as a result 
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of weight economies achieved in the 
system by RCA engineers. For ex- 
ample, the orientation-control system, 
consisting of electrical circuitry and 
250 turns of wire around the outside 
of the satellite, weighs less than 2 Ib. 

In the Tiros II ground system, a 
major change from Tiros I| is the 
transfer of one of the two primary 
ground stations from Hawaii to the 
Pacific Missile Range, with installa- 
tions at Point Mugu, southern Cali- 
fornia, and on San Nicolas Island, 
55 miles offshore. This new station, 
operated by the U.S. Navy for NASA, 
is paired with the second main ground 
station operated by RCA and the 
U. S. Army Signal Corps at Fort 
Monmouth, N. J. Through the two 
stations, Tiros II is programmed for 
television picture-taking and infrared 
measurements during the next orbit, 
and is commanded to “read out” the 
data when it comes again within 
range of one of the stations. This 
communications range is roughly 1500 
miles. The main ground system is 
backed up with an additional ground 
station at RCA’s Space Center at 
Princeton. 

The remainder of the Tiros II 
system includes those elements which 
proved so successful in Tiros I. 
Among them are two specially-de- 
veloped RCA television cameras—a 
wide-angle }-in. Vidicon camera which 
photographs an area of about 700 
miles on each side, and a narrow-angle 
version of the same camera which 
views an area about 80 miles on each 
side. 

Associated with the cameras are 
two specially-designed RCA television 
tape recorders to store the TV pic- 
tures until the ‘‘readout’’ command 
is received from a ground station. 
Each of the tape recorders stores up 
to 32 remote TV pictures from its 
associated camera on each orbit. 
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Like Tiros I, Tiros II is expected to 
operate for about three months, 
after which period its cameras will be 
pointing earthward in darkness. Dur- 
ing this operating period, the satellite 
will complete approximately 1,300 
orbits. 

Among the other major features 
common to both the Tiros II and 
Tiros I satellites are electronic clocks 
to control the timing of operations 
by the cameras, tape recorders, and 
infrared systems during each orbit; 
a coating of more than 9,000 solar 
cells on the top and sides of the satel- 
lite to convert solar energy to elec- 
trical energy for operation of the 
electronic systems; “‘yo-yo” weights 
which will be used as the satellite 
enters orbit to slow its spin rate from 
120 revolutions per minute down to 
12 rpm. rotation desired for operations 
in orbit; and five pairs of solid fuel 
spin-up rockets to restore spin mo- 
mentum as required during operations. 


Cryogenic Memory Plane.—Suc- 
cessful fabrication of a cryogenic thin 
film memory plane the size of a 
postage stamp and the development 
of automatic control techniques ca- 
pable of duplicating it many times 
have been announced by the Interna- 
tional Business Machines Corporation. 

This is the first time that scientists 
have used automatic control tech- 
niques to duplicate such a cryogenic 
unit. The memory plane, consisting 
of 135-cryotron devices built up in 
19 layers of material, was fabricated 
by scientists at IBM’s Federal Sys- 
tems Division laboratory in Kingston, 
N. Y., working in conjunction with 
IBM Research. 

In the device reported by IBM, 
three cryotrons form a memory cell 
which combines storage with an 
elementary logic operation. Because 
of this combination of functions, 
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sophisticated operations not practical 
before may be realized. For example, 
scientists expect this concept may be 
used in an “associative memory,” 
which simultaneously searches all 
memory compartments and thereby 
speeds access to any stored informa- 
tion. 

Key to this major advance in 
cryogenic device development is the 
special technique which permits ac- 
curate duplication of devices. By 
means of this technique, microscop- 
ically thin layers of metals and 
insulating materials are automatically 
deposited on a glass substrate. 

The equipment used for deposition 
allows each layer of a metal or insula- 
tor to besequentially deposited through 
17 microscopically adjusted masks, 
or perforated metal sheets. The 
masks are changed automatically like 
records in a juke box and are held 
in a large metal cylinder operating 
under high vacuum. Once the masks 
have been properly aligned, the process 
automatically produces duplicate su- 
perconducting memory planes with 
similar electrical and mechanical 
characteristics. 

Forty separate pieces or “‘bits” 
of information are stored in 120 
cryotrons in the new memory plane. 
Of the remaining 15 cryotrons, 10 
permit access to the stored bits of 
information; the other five are ‘‘in- 
line” cryotrons which switch bits 
of information from one memory 
plane to another. These newer “‘in- 
line’ cryotrons operate faster than 
conventional cryotrons in communi- 
cating between the memory planes. 
They also permit faster communica- 
tion between the memory planes and 
the arithmetic portion of the com- 
puter. Research is in progress to 
refine the thin-film fabrication tech- 
niques to produce cryogenic devices 
of even greater speeds. 
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Interchange of Computer Programs 
Using Plain English Words.—The 
Radio Corporation of America and 
Remington Rand have announced 
that for the first time computer 
programs utilizing plain English have 
been successfully interchanged be- 
tween data processing systems of 
different manufacturers. 

The report followed a_ two-day 
demonstration by the Remington 
Rand UNIVAC Division of the Sperry 
Rand Corporation and the RCA 
Electronic Data Processing Division 
before the Executive Committee of 
the Conference on Data Systems 
Languages (CODASYL), a govern- 
ment-sponsored group formed to 


achieve computer compatibility. The 
demonstration was an important mile- 
stone in the progress of data proc- 
essing, since the practicability of a 
Common Business Oriented Language 
(COBOL) as an effective program- 
ming language for business-type prob- 


lems was established beyond question. 

In these demonstrations two operat- 
ing problems, one from a government 
agency and one from a large indus- 
trial company, were written by RCA 
and Remington Rand utilizing COBOL. 
Each program was automatically fed 
into the data processing equipment 
of the manufacturer that prepared the 
program. Then, with minor changes, 
it was automatically prepared for the 
other manufacturer’s equipment. In 
both cases the results were outstanding. 

The necessity for preparing costly 
and time-consuming programs, or 
detailed instructions, for each make 
of computer has been a major problem 
facing the electronic data processing 
industry. 

The urgency of the problem was 
recognized as early as May, 1959, 
when CODASYL was formed under. 
the guidance of the Department of 
Defense. The committee consists of 


a 
ae 
: 
{ 
‘ 
4 
ou 
gat 
: 
} 
q 
q 
i 
: 
: 


Feb., 1961.] 


government and industry represen- 
tatives. 

Spurred by the Defense Depart- 
ment’s leadership, manufacturers went 
to work to develop a common language 
for computers. COBOL is the result 
of this research. 

In the two-day demonstration, Rem- 
ington Rand used the UNIVAC II 
at its Philadelphia Engineering Center. 
On the second day the conference 
moved to the RCA 501 Systems 
Center at Cherry Hill, New Jersey, 
where instructions were fed into the 
RCA computer system. 

The programs were interchanged by 
the two companies and were fed into 
the computers with slight modifica- 
tions. The results achieved were 
the same whether a UNIVAC program 
was being handled by the RCA system 
or vice-versa. 

COBOL basically is a program- 
ming system that employs simple 
English words to instruct a computer 
—as opposed to the complicated 
machine code understood only by 
computer specialists. 

Not only does COBOL serve to 
achieve compatibility between com- 
puters but the programming time 
can be cut in half. 

A compiling system operates in this 
fashion: The computer is program- 
med to accept nouns chosen by the 
user that are related to the business 
involved—for example “‘payroll file,” 
“wage,” “employee number,” or “tax.” 
The computer also is programmed to 
accept some twenty verbs that serve 
as commands—‘‘write,” ‘“‘read,’’ ‘‘di- 
vide,” “‘add,”’ “display,’’ or ‘‘move,”’ 
along with conjunctions that link 
the verbs to the plain English nouns. 
The word “‘if” can tell the computer 
to make comparisons of data and 
execute logical decisions as a result 
of the comparisons. 

The English words are transferred 
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to magnetic tape and fed into the 
computer for action. 

Under the COBOL System, a 
computer order can be written “SUB- 
TRACT DEDUCTIONS FROM 
GROSS GIVING NET.” Under the 
strictly numerical programming sys- 
tem, the same instruction may be fed 
into the computer in this form— 
72-010237-00-600000 and a pair of 
other numbers equally as long in the 
case of the RCA 501 or BOOS49 
S00623 CO0942 for the UNIVAC IT. 


Satellite Determination System.— 
A new method to determine the exact 
location, range, and the number of 
satellites orbiting the Earth, was out- 
lined in a speech by one of the nation’s 
foremost astronautics experts, Lt. 
Robert M. L. Baker Jr., at the Seventh 
Annual Science and Engineering Sym- 
posium, in a speech entitled, “‘Novel 
Orbit-Determination Techniques.”’ 

“It would be of military advan- 
tage,’’ Baker said, “‘to include a self- 
contained, friend-or-foe satellite iden- 
tification device on board a Satellite 
Inspector (SAINT)-type vehicle—a 
device that would not be dependent 
upon a continuous communication 
link with either friendly stations 
on Earth or other satellite stations. 

“Such a device,’” Baker continued, 
“might be based upon the unique 
‘license plate’ carried around by every 
space vehicle—that is, its orbital 
elements. 

“The determination of | satellite 
orbits from satellite observations,” 
Baker added, ‘‘offers significant ad- 
vantages over terrestrial orbit deter- 
mination.” 

Baker outlined these advantages: 

First there would be no ionospheric 
refraction problems. 

Second, radar frequencies could 
be chosen that would greatly handi- 
cap a terrestrial station from moni- 
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toring the active orbit-determination 
procedure. 

Third, knowledge could be main- 
tained of launches and space vehicle 
maneuvers anywhere on Earth. 

“From this point of view,’’ Baker 
concluded, ‘‘the data collected by 
these ‘space reconnaissance’ or ‘in- 
spection’ satellites could be utilized 
directly to augment the capability 
of the National Space Surveillance 
System as well as provide our satellite 
systems with an immediate reaction 
capability.” 


Lens-Studded Satellites Aid World- 
Wide Communications.—The U. S. 
Air Force may depend upon lens- 
studded satellites, now under develop- 
ment at the Air Research and Devel- 
opment Command’s Wright Air De- 
velopment Division, Wright-Patterson 
AFB, Dayton, Ohio as relay links in a 
world-wide communications network. 

The satellite is covered with radio 
frequency lenses mounted over spher- 
ical reflectors. It is designed so that 
no matter which position the satellite 
assumes with regard to the Earth, 
strong reflection of signals will take 
place between ground stations 6000 
miles apart. 

The  lens-reflector arrangement 
achieves a high degree of reflectivity 
with a relatively small satellite. The 
lenses focus the radio frequency 
energy at a point on the reflector, and 
again focus the reflected energy into 
a beam encompassing the Earth. 

This offers the advantage over 
Echo-type satellites of directing all 
reflected radio energy back toward the 
Earth. Without a system to steer the 
signals earthward, much of the energy 
would be scattered into space. 

A lens-reflector satellite the same 
size as the current ‘‘satelloon’”’ (100 ft. 
in diameter) will reflect signals 100 
times more powerfully than does 
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Echo, according to WADD scientists. 
They also theorize that spinning or 
tumbling of the satellite will not 
impair its efficiency. 

First test model of the satellite 
will be about 10 ft. in diameter and 
will be suspended from a high-altitude 
balloon for its initial tests. 

Development of this satellite is part 
of the Air Force’s passive communica- 
tions satellite program, which is being 
directed by ARDC’s Air Force Com- 
mand and Control Development Di- 
vision (L. G. Hanscom Field, Bedford, 
Mass.). 

Passive satellites offer many ad- 
vantages over communications satel- 
lites equipped with electronic repeat- 
ing devices. Useful lifetime of passive 
satellites is nearly indefinite, depend- 
ing only upon the orbit attained, and 
passive vehicles also have a much 
higher degree of reliability since 
the communications electronics are 
all ground-based and are readily 
accessible for maintenance. 

The passive satellites also have 
relatively low development and pro- 
duction costs, can handle large band- 
widths, and reach operational availa- 
bility rapidly. Further, electronic 
countermeasures (jamming) must be 
directed against ground equipment, 
rather than satellite-carried instru- 
ments. In ground-based communica- 
tions equipment, size, weight, and 
power requirements are not critical. 

Unlike the Echo satelloon, the 
WADD satellite will be made rigid 
through a system of chemically- 
inflated ribs. This is necessary to 
prevent deformation of the satellite 
and a resulting unplanned scattering 
of reflected radio energy. 


Explosive Gas Detector.—A new 
instrument for the detection of dan- 
gerous or explosive vapor mixtures 
has been introduced by the Instru- 
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ments and Systems Section of the 
Industrial Equipment Division of 
Engelhard Industries, Inc., of New- 
ark, N. J. Called “the Bilge Bug,” 
the device is primarily intended for 
use aboard boats but may also be used 
advantageously in industrial or other 
locations to detect concentrations of 
gasoline, natural gas and many other 
combustible vapors. 

Consisting of a control and indicator 
panel, a sensor unit and a connecting 
cable assembly, the Bilge Bug is 
supplied for two operating voltages 
of 6 and 12 volts d.c., and draws a 
current of only 1.2 amp. The sensor 
unit is installed in any location where 
combustible mixtures may collect and 
the control unit is mounted in some 
convenient location such as the helm 
of a boat. 

In operation, the switch is turned 
on, and a pilot light indicates proper 
operation and illuminates the face 
of the meter. Within seconds, the 


two platinum wire sensing elements 
will react to the presence of vapors 
and indicate the condition on the 


meter. Following the check, if the 
reading is ‘‘safe,’’ the instrument may 
be turned off and the engine started. 


Self-Extinguishing Epoxies.—De- 
sign and production of two self- 
extinguishing epoxy resins are an- 
nounced by The Dow Chemical 
Company of Midland, Mich. Both 
resins are made self-extinguishing 
through halogen—bromine—substi- 
tuted on the epoxy molecule. They 
are cured or hardened with the common 
epoxy curing agents. 

According to Dow, self-extinguish- 


ing epoxy resins should find use in. 


laminated aircraft structures, glass 
and paper-based laminated electrical 
circuits, filled castings, and in aircraft 
and missile adhesives. Marine paints, 
architectural coatings and _ sheath 
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coatings (for solid or foamed plastics) 
are other likely uses. 

Experimental Resin X-3442 is a 
semi-solid containing about 49 per 
cent bromine and has utility in blends 
with conventional liquid epoxy resins. 
Experimental Resin X-3441.1 is a 
solid resin containing about 19 per 
cent bromine and can be used either 
alone or in blends with conventional 
solid resins. A bromine content of 
15-25 per cent is said to make even 
clear castings self-extinguishing. 

Physical and electrical properties 
of castings and laminates made with 
the new resins are comparable to 
properties developed by conventional 
epoxy resins. Modulus properties 
are reduced at 300°F., but blends of 
self-extinguishing epoxy with Dow’s 
epoxy novolac resin are expected by 
the company to result in higher 
modulus at elevated temperatures. 

Prepreg laminates made with either 
of the self-extinguishing resins meet 
the strict NEMA G-11 physical and 
electrical properties except for 300°F. 
flexural strength. The resins are 
therefore suitable for use in many 
commercial applications. 


“Talking” Through Touch.—A re- 
search psychologist has revealed the 
feasibility of giving pilots another 
channel of communications, for in 
these days of supersonic speeds, eyes 
and ears cannot cope with the mass of 
communications fed into a pilot’s 
brain. Why not, research scientists 
reasoned, give pilots a third channel 
of communications? 

Under the direction of Psychologist 
Richard A. Newman, the Human 
Engineering Laboratory of the Rome 
Air Development Center (RADC) 
undertook the research program to 
add another communication channel 
to the brain. RADC is a sub-unit 
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of the Air Research and Development 
Command. 

The program was revealed in a 
paper prepared by Newman for de- 
livery at ARDC’s Seventh Annual 
Science and Engineering Symposium. 
' The new method involves change 
of vibrations, or frequencies of the 
voice, into mechanical vibrations which 
a pilot can feel by placing a small plate 
in contact with the skin. 

Actually, the plate could be placed 
anywhere on the body. But research 
revealed that an area on the pilot’s 
chest is more receptive to pulses than 
other parts of the body. 

Newman emphasized that the proj- 
ect still is strictly in the research stage. 
He added, however, he felt scientists 
were making substantial progress with 
the technique. 

“So far as application for the Air 
Force is concerned,’’ Newman said, 
“we are thinking of developing the 
touch system for channeling informa- 
tion to a pilot when his hearing or 
sight has been impaired by injury, 
weather, or environmental conditions.” 

Scientists envision numerous con- 
ditions under which the technique 
might prove practical. 

On takeoff, for instance, jet pilots 
perform many manual functions while 
visually checking numerous instru- 
ments and talking with the control. 
tower operator. Use of a _ small 
vibrator on the body would give them 
another means of receiving certain 
types of information. 
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The psychologist pointed out that 
the system could be used as a back-up 
channel when jet noise or other 
technical difficulties disrupt normal 
communications. 

Another situation when this means 
of communication might prove prac- 
tical would be individuals responsible 
for missile firings. As the moment of 
launch approaches, numerous _indi- 
viduals must perform multiple tasks 
almost simultaneously. A third com- 
munications channel would enable 
them to perform more efficiently. 

A third use of the “‘sensor’’ might 
be when an individual would want 
to communicate with another without 
the knowledge of persons nearby. 
A small vibrator attached to the body 
could easily go unnoticed. 

“As in other Air Force research 
projects, this system may have far- 
reaching civilian applications,’’ New- 
man said. 

In the future, touch communica- 
tions might be used to assist deaf 
people to “hear” through their skin. 
This might be done by training them 
to feel the vibrations from voices by 
using a small vibrator. 

Newman says the technique offers 
much promise to individuals asso- 
ciated with the aerospace age. He 
added, however, with typical scien- 
tific caution, that “a more extensive 
research program will be required 
before a fully operable technique 
can be developed.” 
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THE FRANKLIN INSTITUTE 
LABORATORIES FOR RESEARCH AND DEVELOPMENT 


announce opportunities for: 


ELECTRICAL ENGINEERS 
SYSTEMS ENGINEERS 
ORGANIC CHEMISTS 

OPERATIONS RESEARCH SPECIALISTS 
MECHANICAL ENGINEERS 
EXPERIMENTAL PSYCHOLOGISTS 
METALLURGISTS 
INDUSTRIAL ENGINEERS 
PHYSICISTS 
APPLIED MATHEMATICIANS 


to work on challenging problems in research 


Send complete resume to: 
Mr. John E. Christ, Director of Personnel 


THE FRANKLIN INSTITUTE 
Philadelphia 3, Pa. 
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Now in the fourth year of operation, The Franklin Institute Com- 
puting Center has shared its technical know-how, modern electronic 
equipment and experienced personnel with hundreds of progressive 
industries and government agencies across the nation. 


We offer the services of creative people, highly skilled in their respec- 
tive fields and ably trained in the application of these skills to the ever 
expanding area of electronic computers and data processing systems. 


This staff is now available for analysis, system design, programming 
or coding of projects of unlimited scope or context. Input to our large 
scale computer and peripheral equipment is acceptable in any form. 
Results are provided on cards, plastic or metallic tape, and in com- 
pletely edited printed copy. Our extensive library of automatic coding, 
engineering, data processing and mathematical routines is available to 
all users, and machine time is provided with or without the services 
of programming personnel. 


THE FRANKLIN INSTITUTE 
Computing Center 


20th Street and Benjamin Franklin Parkway 
Philadelphia 3, Pa. LOcust 4-3600 
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